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ABSTRACT 


The  overall  purpose  of  this  work  is  to  develop  a sub- 
sonic, turbulent,  diffuse  electric  discharge  operating  at 
near  atmospheric  pressure  for  use  as  the  pumping  mechanism 
in  an  Electric  Discharge  Convection  Laser  (EDCL) . The 
work  was  done  in  two  phases.  The  first  consisted  of  a 
data  collection  phase  during  which  the  variation  of  dis- 
charge parameters  was  studied  to  determine  their  effects 
on  the  current-voltage  characteristics  of  the  multi-pin 
electrode  arrangements  used.  The  parameters  studied 
included  gap  length,  flow  speed,  turbulence  intensity, 
electrode  orientation,  rate  of  voltage  increase,  and  gas 
density  through  the  addition  of  a diffuser.  The  second 
phase  consisted  of  developing  an  interactive  model  which 
correctly  reflects  the  dependences  observed  in  phase  one. 

A space-charge  model  resulted  which  correctly  explains 

most  major  observed  trends  in  gross  discharge  properties. 

* 

Discharge  stability  is  also  discussed  in  terms  of  the  model. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION - 12 

A.  GENERAL 12 

B.  THE  MODEL 14 

II.  EXPERIMENTAL  APPARATUS 19 

III.  EXPERIMENTAL  OBJECTIVES  AND  MODELING 

APPROACH - 23 

A.  GENERAL - 23 

B.  MODELING  APPROACH  - 24 

C.  PREVIEW  OF  EXPERIMENTAL/MODEL 

CORRELATIONS  26 

IV.  SPACE  CHARGE  MODEL  28 

A.  BREAKDOWN  VOLTAGE  (Vfo)  29 

1.  Uniform  Field  (parallel  plate)  29 

2.  Non-Uniform  Fields  33 

a.  No-Flow 33 

(1)  Spatial  Development  of 

Virtual  Anode  37 

b.  Gas  Convection/Flow  Orientation 

Effects  Analysis  41 

(1)  Space-Charge-Convection- 

Effects  Model 42 

(2)  Experimental  Inputs  and 

Basic  Model  45 

(3)  Maximum  Voltage  Ratio  Due 

to  Convection 47 


5 


(4)  Minimum  Voltage  Ratio 49 

(5)  Turbulence  Effects  on 

Breakdown  Voltage  56 

(6)  Cross-Flow  Stability  and 

Field  Interactions  57 

(7)  Effects  of  Charged  Particle 
Density  on  Recombination 

Times  60 

(8)  Reverse  Flow  Convection 

Effects  67 

(9)  Analytic  Considerations  in 

I = f (U)  — 70 

B.  GASDYNAMIC  EFFECTS  ON  DISCHARGE  CURRENT 

(NON-UNIFORM  FIELDS)  72 

1.  Pure  Convection  Effects  on  Current  76 

2.  Current-Voltage  Characteristics  77 

3.  Diffuser/Density  Effects  82 

4.  Turbulence  Effects  84 

5.  Maximum  Current/Discharge  Stability  90 

V.  CONCLUSIONS  93 

VI.  RECOMMENDATIONS  - 95 

ILLUSTRATIONS  - 98 

APPENDIX  A VARIABLE  VOLTAGE  RATE  RESULTS  129 

APPENDIX  B EMPIRICAL  FIT  TO  PUBLISHED  BREAKDOWN 

DATA 132 

BIBLIOGRAPHY  — - 143 

INITIAL  DISTRIBUTION  LIST  - — 145 


6 


LIST  OF  ILLUSTRATIONS 


F igure 

1.  Sample  Electric  Field  Distributions  98 

2.  Electrode  Design  99 

3.  Turbulence  Generating  Grids  100 

4.  Power  Supplies 101 

5.  Variable  Speed  Transmission  102 

6.  Parallel  Plate  Space  Charge  Effects  103 

7.  Comparison  of  Space  Charge  Model  to  Data 

for  E^ 104 

8.  Breakdown  Voltage  Density  Dependence  105 

9.  Schematic  of  Flow  Orientations  106 

10.  Breakdown  Voltage  Data  in  Ordinary  Flow  Mode  --  107 

11.  Breakdown  Voltage  Data  in  Reverse/Cross- 

Flow  Modes  108 

12.  Point-Plane  Field  Distribution  with  Space 

Charge 109 

13.  Breakdown  Voltage  vs.  Gap  Length 110 

14.  Predicted /Ob served  V^/V^n)may  vs.  Gap  Length  --  111 

15.  Slip  Parameter  at  E = 20  kV/cm 112 

16.  Slip  Parameter  at  Anode  Surface,  E = E 113 

P 

17.  Development  and  "Blowout"  of  Anode  Extension  --  114 

18.  Schematic  of  Observed  Pin  Tip  Glow  Patterns  115 

19.  Photographs  of  the  Discharge 116 


F igure 

20.  Photographs  of  the  Discharge  117 

21.  Schematic  of  Cross-Flow 118 

22.  Recombination  Lengths  in  Cross-Flow  Mode  119 

23.  vs.  Velocity  as  Predicted  by  Model 120 

24.  Reverse  Flow  Schematic  121 

25.  Avalanche  Travel  for  Various  Conditions  122 

26.  The  Current  Equation 123 

27.  I-V  Characteristic  Velocity  Shift  124 

28.  Vb  vs.  Velocity  Summary  Plot 125 

29.  Ionization  Efficiency  of  Air 126 

30.  I-V  Characteristic  Density  Shift  127 

31.  Breakdown  Current  vs.  Mass  Flow  Rate 128 

A-l  Variable  Voltage  Rate  Data 131 

B-l  Breakdown  Power  vs.  Density  Ratio  135 

B-2  Breakdown  Voltage  vs.  Density  Ratio  136 

B-3  Breakdown  Current  vs.  Density  Ratio  137 

B-4  Multi-Pin  Specific  Power  vs.  Gap  Length  Data  --  138 

B-5  V^  vs.  Velocity  Data  with  Screen  # 9 139 

B-6  1^  vs.  Velocity  Data  with  Screen  # 9 140 

B-7  P^/cm  vs.  Velocity  Curves  with  Screen  # 9 141 

B-8  I.  vs.  Gap  Length  at  Various  Velocities/ 

Screen  # 9 142 


8 


LIST  OF  SYMBOLS 


Symbol 

A Stoletow  Constant 

a Constant 

B Stoletow  Constant 

C Constant 

Cl  Constant 

C2  Constant 

c Constant 

d Gap  Length  (actual) 

d'  Effective  Gap  Length 

E Electric  Field  Strength 

E'  Effective  Electric  Field  Strength 

E^’  Ionization  Field  Strength 

e Electronic  Charge 

f Frequency  of  Current  Pulses 

f(  ) Function  of  ( ) 

I Current  (total) 

I Breakdown  Current 

D 

I No -Flow  Current 

o 

i Avalanche  Current 

j Current  Density 


J 


9 


k 


| 


M 

N 

n 

o 

P 

P 

P 

pi 

R 

S 

t 

U 

V 

v 


vd 


V 


b 


V 

c 

X 


x 

X 

rec 


Charged  Particle  Mobility 
Corona  Initiation  Voltage  = V£ 

Number  of  Avalanches 
Particle  Density 

Subscript  for  No-Flow  or  STP  Value 
Pressure  (total) 

Amplitude  of  Fluctuating  Pressure  Component 
Power 

Breakdown  Power 

Tip  Radius  of  Pin  Anode 

Slip  Parameter 

Time 

Flow  Velocity 

Applied  Voltage 

Total  Velocity 

Drift  Velocity 

Breakdown  Voltage 

Corona  Initiation  Voltage 

Anode  Extension  due  to  Convection  (Ut  ) 

Axial  Distance  from  Anode  Surface 

Recombination  Distance 

Ionization  Distance 

Ionization  Coefficient 


. 


10 


a 

r 


Recombination  Coefficient 


Rate  of  Voltage  Increase 

Anode  Extension 

No -Flow  Anode  Extension 

Mean  Free  Path 

Gas  Density 

Standard  Gas  Density 

Turbulent  Frequency 

Summation 

Integral 

Ionization  Efficiency  ( ai/E) 
Characteristic  Time  Interval 
Proportional  to/Approximately  Equal  to 
Order  of  Magnitude  of 


I.  INTRODUCTION 


A.  GENERAL 

This  work  is  the  continuation  of  an  ongoing  project  at 
the  Naval  Postgraduate  School  (NPS)  where  the  feasibility 
of  constructing  a high-energy,  high  pressure,  continuous- 
wave  Electric  Discharge  Convection  Laser  (EDCL)  is  being 
investigated. 

Laser  systems  of  this  type  currently  under  investigation 
are  limited  in  power  output  by  the  amount  of  electrical 
power  which  can  be  coupled  into  the  laser  (Ref.  1).  Since 
laser  gain  is  directly  proportional  to  the  population 
inversion  created  by  electrical  pumping,  the  importance  of 
optimum  conditions  in  the  discharge  region  can  readily  be 
seen. 

Previous  work  on  this  project  at  NPS  vividly  demonstrated 
the  importance  of  flow  conditioning  on  discharge  performance. 
Power  increases  in  excess  of  250  times  the  normal  no-flow 
conditions  were  demonstrated  by  flowing  air  rapidly  through 
the  discharge  region  with  proper  conditioning  by  turbulence 
generating  screens  (Ref.  2). 

A brief  general  look  at  the  physical  situation  pre- 
vailing in  the  discharge  section  will  be  presented  here 
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to  review  the  state  of  the  art  and  to  aid  in  interpretation 
of  observed  electrical  characteristics.  The  basic  geometric 
design  of  the  electrodes  is  a common  one  often  used  in 
electric  discharge  work,  the  positive  point-to-plane . This 
design  results  in  an  electric  field  distribution  in  the 
absence  of  space  charge  which  is  proportional  to  1/x,  where 
x is  the  distance  from  the  point.  This  inverse  field 
dependence  is  common  to  many  non-uniform  field  electrode 
geometries  including  the  hyperboloidal  point-to-plane, 
concentric  cylinders,  wire-to-plane,  sphere-to-plane,  and 
many  others,  so  the  analysis  presented  here  can  be  applied 
to  a number  of  technically  interesting  cases. 

The  field  inverse  dependence  on  distance  from  an  elec- 
trode results  in  the  division  of  the  electrode  gap  into 
distinct  regions  of  very  high  field  where  x is  small,  and 
very  low  field  when  x is  relatively  large.  The  smaller 
the  radius  of  the  electrode  involved,  the  more  pronounced 
the  division.  Figure  (1)  illustrates  this  effect  for  two 
geometries . 

As  a result  of  the  above,  as  voltage  is  increased  across 
the  gap,  ionization  will  begin  in  the  high  field  region, 
confining  the  discharge  activity  to  the  immediate  vicinity 
of  the  pin  where  free  electrons  are  able  to  gain  enough 
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energy  to  ionize  by  collision.  Just  prior  to  complete 
breakdown  of  the  gap,  there  exists  a very  small  region  of 
intense  ionization  which  appears  as  a bright  glow  at  the 
pin  tip  (in  a point- to -plane  discharge),  and  a space  charge 
region  of  larger  proportions,  caused  by  less  intense  ioniza- 
tion and  electrostatic  repulsion  of  the  positive  ions  formed 
in  the  ionization  region  by  the  positive  pin.  Beyond  this 
space  charge  region  lies  a quasi-neutral  region  which  con- 
sists of  most  of  the  low  field  region  of  Figure  1.  The 
cathode,  while  lying  well  outside  the  high  field  region, 
does  not  play  an  important  role  in  the  dynamic  processes 
taking  place  near  the  anode,  and  so  will  not  be  discussed 
in  this  work,  except  in  its  role  at  final  breakdown. 

B . THE  MODEL 

The  analysis  which  follows  is  based  on  the  ability 
of  the  flow  field  to  alter  the  shape  and  extent  of  the 
space  charge  region  by  convection  and  turbulence.  The 
problem  is  complex  and  involves  many  variables  including: 

(1)  tip  radius  (R) 

(2)  gap  length  (d) 

(3)  mean  gas  density  ( p)  and  composition 

(4)  mean  flow  velocity  (U) 
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(5)  electric/flow  field  coupling  mechanisms  (E,U,k) 

(6)  time  dependent  fluctuations  of  3 and  4 above 
due  to  the  introduction  of  turbulence,  U(t), 


p(t) 

(7)  ionization  and  charged  particle  dynamics  (a  , 
a , X , e,  n) 

The  result  is  that  the  easily  measurable  gross  discharge 
electrical  properties  such  as  voltage  and  current  are 
functions  of  these  variables,  i.e.; 


I~f 1(R»  d,p  , U,  E,  k,  t) 
V=f 2 (R , d,p  , U,  E,  k,  t) 


(1) 


There  are  many  characteristic  times  associated  with  the 
variables  cited,  such  as  gas  transit  times,  characteristic 
times  associated  with  turbulent  fluctuations  of  density 
and  velocity,  characteristic  ionization  and  recombination 
times,  streamer  propagation  and  growth  times,  and  many 
others  (Table  I).  The  times  cited  above,  however,  will  be 
of  primary  interest  in  the  discussion  which  follows. 
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TABLE  I 


Process 

Ionization 

Excitation 


Characteristic  Time  (Sec)  Reference 


Energy  Transfer  from 
Electrons  to  Atoms  — 


Maxwellian  Distrib. 
of  Atoms  


Recombination  

Current  Pulsing  

Streamer  Propagation 
Turbulence  Freq  


Rate  of  Voltage  Increase 
( V/  y ) 


iq-12  

------  in 

iu  — 

io-12 

in 

xu  — 

1 0-5  _ 1 0-6  

in 

XU  XU  — — — - 

10-8  10-10  

in 

XU  “ XU  — — — - 

10-5  10-6  

XU  XU  — — — — 

10— 4 i0-5  

. A 

XU  XU 

10*®  - -- - 

------  A 

10"1-10-5  

3 

- 0-2.5(10  ) 

(NOTE:  The  times  from  Ref.  10  are  for  an  arc  discharge.) 


(present 

work) 


Since  the  model  presented  in  the  following  sections 
depends  on  the  ability  of  the  flow  field  to  couple  with  the 
charged  particles  produced  in  the  discharge  region,  it  is 
important  to  understand  how  this  coupling  may  take  place. 
The  three  types  of  charged  particles  of  primary  interest 
are,  of  course,  positive  and  negative  ions,  and  free  elec- 
trons, The  bulk  of  the  gas  remains  un-ionized  over  the 
major  portion  of  the  discharge. 
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Charged  particles  in  the  gas  under  the  influence  of 
an  electric  field  acquire  a drift  velocity  (V^=kE)  which 
depends  on  the  gas  composition,  pressure,  and  electric 
field  strength.  If  the  gas  is  flowing,  then  the  particle 
velocity  (in  vector  form)  relative  to  the  stationary 
electrodes  is: 

vtotal-  kE  + 0 (2) 

The  sign  convention  used  here  will  be  that  E,  U and  x are 
positive  from  anode  to  cathode.  The  k in  Eq.  (2)  is  the 
charged  particle  mobility,  which  can  be  expressed  in  terms 
of  the  parameters  listed  as  k i e X/mv,  where  m is  the 
particle  mass,  X (the  particle  mean  free  path)  is  inversely 
proportional  to  gas  pressure,  and  v is  the  random  thermal 
velocity  of  the  particle.  Equation  (2)  can  be  rewritten 
as  vtota^  = kE(l  + U/kE) , where  the  quantity  U/kE  will  be 
referred  to  as  the  Mslip  parameter,  S."  Electrons,  due  to 
their  small  mass,  have  a very  high  mobility,  and  generally, 
the  slip  parameter  for  electrons  is  very  small,  and  gas 
convection  is  not  significant;  but  ions  are  much  more 
massive,  so  k is  low  and  ion  drift  velocities  can  be  of 
the  same  order  of  magnitude  as  U,  even  in  the  high  field 
regions  of  the  discharge,  resulting  in  significant  ion 


1 


17 


I 


A 


It 


relocation  by  convection  within  the  characteristic  times 
mentioned  above.  Since  high  densities  of  charged  particles 
produce  electrostatic  fields  of  their  own,  which  can  easily 
reach  the  same  order  of  magnitude  as  the  applied  field,  the 
electrons  are  affected  by  convection  indirectly  through 
the  alteration  of  the  electric  field  structure  by  the  move- 
ment of  high  density  ionic  space  charges. 

Visual  observations  of  a discharge  vividly  demonstrate 
the  ability  of  highly  turbulent  flow  to  affect  the  space 
charge.  With  high  velocity  and  no  turbulence,  the  glow 
extending  from  the  tips  of  the  pins  toward  the  cathode  is 
very  bright  and  steady,  while,  when  high  turbulence  is 
introduced,  the  bright  portion  of  the  glow  can  be  seen  to 
radically  "whip"  around,  indicating  the  high  degree  of 
influence  that  the  fluctuating  velocity  components  have 
on  the  charges  present  in  the  flow. 
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II.  EXPERIMENTAL  APPARATUS 


The  primary  test  apparatus  used  in  this  work  consists 
of  a plexiglass  wind  tunnel  with  a test  section  in  which 
various  electrodes  and  turbulence  generating  grids  were 
mounted.  Electrode  spacing  is  adjustable  and  the  electrode 
orientation  with  respect  to  the  flow  field  could  be  changed, 
permitting  much  flexibility  in  the  examination  of  flow  and 
electric  field  interactions  for  a variety  of  conditions  and 
with  several  levels  of  turbulence.  Later  in  the  testing, 
a new  wind  tunnel  was  added  which  was  equipped  with  a 
diffuser  to  permit  operation  at  higher  flow  velocities  and 
reduced  density  in  the  discharge  section  (Ref.  3) . 

The  electrodes  are  of  a design  common  to  many  electric- 
discharge  lasers  (Fig.  2).  The  anode  consists  of  a pin 
rack  with  three  rows  of  pins,  eleven  to  thirteen  pins  per 
row.  The  pins  were  made  of  steel,  with  horizontal  and 
vertical  separation  of  one  centimeter.  The  cathodes  used 
varied,  because  there  was  a need  to  search  for  the  most 
effective  and  structurally  sound  arrangement.  All  the 
cathodes  used,  however,  were  geometrically  similar,  and 
simulated  a plane  with  minimum  solid  area  to  reduce  flow 
blockage  as  much  as  possible.  Structural  integrity  became 
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an  important  consideration  as  testing  progressed,  since 
flow  velocities  exceeded  180  meters  per  second  and  turbulence 
levels  introduced  were  relatively  high  (Ref.  3). 

All  tests  were  performed  in  an  air  flow  which  was 

supplied  by  a Carrier  three-stage  Centrifugal  Compressor 

3 

with  a maximum  inlet  capacity  of  4000  ft  /min,  and  a maximum 
pressure  ratio  of  two.  Flow  velocity  was  regulated  through 
a series  of  valves,  and  flow  total  temperature  was  kept 
below  100°F  by  passing  the  air  through  an  air  cooler  prior 
to  reaching  the  test  section.  Flow  velocity  was  measured 
with  pitot  or  pitot-static  tubes  which  were  mounted  in  or 
downstream  of  the  test  section,  depending  on  the  tunnel 
being  used.  Turbulence  was  introduced  by  mounting  one  or 
more  grids  upstream  of  the  discharge  section.  Some  grids 
which  were  used  in  previous  work  were  retested  for  verifi- 
cation and  validation  of  systems  operation  (Ref.  3) , and 
several  new  designs  were  tested  (Fig.  3). 

Turbulence  spectra  were  measured  with  a Thermo  Systems, 

Inc.  hot-wire  anemometer,  which  consisted  of  a 1051-2  monitor 
and  power  supply,  1054A-30  anemometer  module,  and  locally 
manufactured  hot-wire  probes.  A detailed  account  of  turbu- 
lence measurement  equipment  and  techniques  can  be  found  in 
Reference  (3).  4 

J 


20 


Two  power  supplies  were  used  to  power  the  electric  dis- 
charge. The  first  was  a Sorensen  Model  1030-20  High  Voltage 


D.C.  power  supply,  which  is  rated  at  30  kV  D.C.  and  20  ma 
with  a 27o  ripple.  It  was  soon  discovered  that  the  maximum 
gap  length  was  severely  limited  by  this  supply,  so  a Sorensen 
BETA  Model  High  Voltage  Supply  was  obtained  which  was  rated 
at  60  kV  D.C.  and  50  ma  (Fig.  4).  With  this  power  supply 
it  was  hoped  to  increase  the  maximum  gap  length  which  could 
undergo  breakdown.  Due  to  uncorrectible  internal  problems, 
the  maximum  voltage  and  current  obtainable  was  34  kV  and  19  ma, 
so  reliable  tests  were  limited  to  less  than  4 cm  gap  length. 
Both  power  supplies  were  equipped  with  trip  controls  which 
work  on  either  voltage  or  current.  The  voltage  was  initially 
measured  by  a Sensitive  Research  electrostatic  voltmeter 
with  a maximum  range  of  40  kV  and  an  internal  impedance  of 
5.10^  ohms.  The  readings  of  the  power  supply  voltmeters 
were  compared  and  were  found  to  be  in  excellent  agreement, 
so  they  were  used  exclusively  throughout  the  remaining  tests. 
The  current  was  measured  with  a variety  of  ammeters,  ranging 
from  microamp  scales  to  a 200  ma  scale.  A meter  panel  was 
set  up  with  three  ammeters  in  parallel,  one  0-1  ma  range, 
one  0-50  ma,  and  one  0-200  ma.  The  meters  were  wired 
through  a multi-position  '*Make-before-break"  switch  so  that 


21 


J 


as  the  range  of  one  meter  was  exceeded,  the  next  one  could 
be  switched  in  without  circuit  interruption.  A Simpson 
volt-ammeter  was  also  wired  in  series  with  the  ammeter 


panel  for  use  in  the  microampere  range  and  as  a check  against 
the  readings  of  the  panel  meters.  This  redundancy  permitted 
immediate  detection  of  erroneous  current  readings,  and  re- 
placement of  defective  meters.  All  wiring  was  done  with 
high  voltage  leads  which  were  fabricated  prior  to  the  initial 
test  phase.  Connections  were  made  through  highly  polished 
brass  spheres  which  were  about  1\  inches  in  diameter.  The 
high  voltage  side  from  the  power  supply  outlet  to  the  anode 
was  immersed  in  a high  dielectric  oil  bath  to  prevent 
voltage  leaks  at  the  ball  connection.  Great  care  was  taken 
to  ensure  that  all  equipment  was  well  grounded  at  all  times, 
through  the  laboratory  grounding  system. 

One  series  of  tests  performed  required  increasing  the 
applied  voltage  across  the  gap  at  various  rates.  A constant- 
speed  24  volt  D.C.  motor,  powered  by  an  adjustable  output 
D.C.  Power  Supply  was  connected  to  the  Sorensen  30  kV 
power  supply  voltage  control  shaft  through  a Metron  Variable- 
Speed  transmission  (Fig.  5).  This  arrangement  permitted  the 
applied  voltage  to  be  increased  across  the  gap  at  rates  from 
0-20,000  volts/second. 

22 


W 


III.  EXPERIMENTAL  OBJECTIVES  AND  MODELING  APPROACH 


A.  GENERAL 

A brief  preview  of  the  experimental  work  and  an  overall 
perspective  of  the  model  development  should  serve  to  help 
keep  sight  of  the  "big  picture"  while  sorting  out  the  de- 
tails of  the  individual  pieces. 

The  experimental  work  has  been  divided  into  two  main 
categories.  First,  tests  were  run  varying  individual  dis- 
charge parameters,  including  gap  length,  flow  velocity,  gas 
density,  turbulence  intensity,  cathode  configuration,  and 
electrode  orientation.  The  purpose  of  these  tests  was  to 
establish  a data  base  for  use  in  analyzing  the  individual 
contributions  of  various  parameters  to  discharge  power, 
and  to  attempt  to  formulate  an  interactive  physical  model 
based  on  the  observed  relations.  The  second  phase  of  test- 
ing consisted  of  observing  the  effects  of  increasing  applied 
voltage  across  the  discharge  gap  at  various  rates  to  deter- 
mine what  effect  this  would  have  on  the  breakdown  properties 
of  the  gap.  The  rate  of  voltage  increase  was  constant  for 
each  run,  but  it  varied  from  run  to  run.  These  tests  also 
added  to  the  data  base  which  was  being  established,  but  the 
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specific  purpose  was  to  determine  if  any  type  of  resonance 
effect  might  be  detected  between  a particular  rate  of  in- 
crease, and  the  turbulence  spectra  which  were  being  intro- 
duced since  the  characteristic  time  scales  of  turbulence 
and  the  rates  of  increase  were  of  the  same  order  of  magnitude. 

B . MODELING  APPROACH 

A space-charge  model  development,  which  begins  in 
Section  IV,  is  based  on  the  dynamic  processes  which  occur 
when  the  gap  is  very  near  its  breakdown  potential.  Under 
these  conditions,  a positive  point- to-plane  gap  in  an 
electronegative  gas  such  as  air,  experiences  a pulsating 
current  due  to  the  formation  of  pre-breakdown  streamers 
whose  magnitudes  vary  in  a statistical  way.  The  frequency 
of  these  pulses  has  been  measured  on  the  order  of  10^  to 
10^  Hz  (Ref.  4) . The  model  proposed  here  looks  at  the 
discharge  from  the  point  of  view  of  an  individual  event, 
such  as  the  propagation  of  a pre-breakdown  streamer,  and 
the  situation  that  prevails  in  the  discharge  at  the  time 
and  location  of  its  occurence.  For  example,  in  a steady, 
non-flowing  point-plane  discharge  near  breakdown  potential, 
the  streamers  which  cyclically  form,  all  "see"  generally 
the  same  "picture"  (field  structure,  space  charge  distribu- 
tion, etc.)  as  the  preceding  ones  did.  Much  activity  takes 
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place  between  these  events,  but  at  each  period  the  situa- 
tion  at  a given  location  is  very  similar.  Now,  if  some- 
thing is  done  to  alter  the  picture  that  the  streamers  see, 
such  as  a reduction  of  density,  or  flowing  gas  through  the 
gap,  or  the  introduction  of  turbulence,  the  "picture"  pre- 
sented to  each  streamer  may  be  very  different  from  the 
steady,  no-flow  case.  The  change  in  the  local  situation 
can  be  the  result  of  two  separate  classes  of  influencing 
factors,  one  steady-state,  such  as  a permanent  change  in 
gas  density;  and  the  other,  a time-dependent  change,  such 
as  with  the  introduction  of  convection  or  turbulence,  where 
the  total  change  seen  by  the  next  streamer  depends  on  how 
much  time  the  effect  has  to  work  on  the  discharge  between 
streamers.  The  change  in  this  picture  over  the  original 
case,  then,  is  what  the  model  assumes  determines  the  changes 
observed  in  gross  discharge  properties.  Of  course,  both 
classes  of  influence  can  be  operative  at  once,  and  then 
their  individual  effects  must  be  sorted  out. 

To  summarize,  the  model  was  developed  by  isolating  a 
particular  influence  factor,  determining  its  effect  on  the 
"picture"  seen  by  successive  streamers,  then  attempting  to 
combine  the  results  of  all  such  influences  in  an  interactive 
way  to  explain  observed  changes  in  gross  discharge 
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characteristics.  The  ultimate  goal  of  this  study  is  to 
develop  a method  of  approach  for  optimizing  the  discharge 
for  use  as  a pumping  mechanism  for  the  EDCL. 


C.  PREVIEW  OF  EXPERIMENTAL/MODEL  CORRELATIONS 

The  first  part  of  Section  IV  introduces  the  basic  space- 
charge  model  and  illustrates  how  a reduction  in  density  by 
addition  of  a diffuser  (Ref.  3)  alters  the  picture  presented 
to  succeeding  streamers.  Successful  development  of  a simple 
relation  for  the  density  dependence  of  gap  breakdown  voltage 
in  a non-uniform  field  evolves  from  fitting  existing  data 
to  an  equation  form  resulting  from  the  space-charge  model. 
This  lends  credibility  to  the  basic  concept. 

The  next  part  shows  an  analysis  of  the  influence  of 
convection  on  the  discharge  with  various  orientations  of 
electric  and  flow  fields.  With  the  aid  of  the  space-charge 
model  again,  the  change  in  "picture"  seen  by  successive 
streamers  is  determined,  and  various  maxima  and  minima  for 
breakdown  voltages  are  calculated  for  variations  of  velocity 
and  orientation.  Again  the  results  obtained  agree  with 
experimentally  observed  values,  lending  further  credibility 
to  the  proposed  model  concept.  Various  instabilities  and 
observed  visual  characteristics  are  also  discussed  in  terms 
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of  the  model,  with  general  qualitative  agreement  between 
model  predictions  and  observations. 

In  the  second  half  of  the  development,  velocity,  density, 
and  turbulence  effects  on  discharge  current  and  the  inter- 
relation between  the  current  and  voltage  are  discussed.  A 
current  equation  evolves  from  the  space-charge  model  which 
properly  accounts  for  most  major  observed  trends  in  dis- 
charge current  due  to  convection,  turbulence,  and  density 
reduction  in  various  flow  orientations.  Some  stability 
conditions  are  mentioned  which  are  particularly  relevant 
to  the  discussion,  but  no  attempt  is  made  to  give  a detailed 
account  of  stability  effects. 

The  four  basic  discharge  conditions  used  for  comparison 
of  the  space-charge  model  predictions  are  the  following: 

(1)  No  diffuser,  no  turbulence 

(2)  No  diffuser,  with  turbulence 

(3)  Diffuser  on,  no  turbulence 

(4)  Diffuser  on,  with  turbulence 

Test  results  appear  to  bear  out  the  predicted  trends  for 
each  of  the  above  cases.  Detailed  design  data  on  the  dis- 
charge apparatus,  including  diffuser  and  turbulence 
generators  are  given  in  Reference  (3) . 
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IV.  SPACE  CHARGE  MODEL 


Numerous  studies  have  been  done  on  various  aspects  of 
flowing  gaseous  discharges  in  order  to  determine  the  effects 
of  convection  and/or  turbulence  on  the  strength  and  stability 
of  such  discharges  (Refs.  2,  3,  5,  6). 

As  pointed  out  in  the  Introduction,  understanding  the 
basic  processes  in  a non-uniform,  flowing  discharge  is 
essential  to  the  improvement  of  the  ultimate  potential  of 
the  Electric  Discharge  Convection  Laser.  Present  discharge 
stability  problems  restrict  practical  operation  of  this 
type  of  laser  to  the  realm  of  100-200  torr  operating  pres- 
sure (Ref.  5).  This  pressure  limitation,  or  more  accurately, 
density  limitation,  in  turn  limits  maximum  power  input  to 
the  gas,  and,  consequently,  restricts  laser  power  output. 

In  the  following  discussion,  an  attempt  is  made  to  shed 
some  light  on  the  complex  problem  of  discharge  stability 
and  power  enhancement  in  a non-uniform,  subsonic,  near- 
atmospheric  flowing  discharge.  Using  some  well-known 
principles  of  gaseous  breakdown  in  uniform  fields  and  a 
few  basic  gas  dynamic  ideas,  it  is  possible  to,  at  least 
qualitatively,  explain  many  diverse  observed  phenomena  of 


the  particular  discharge  sections  used  in  this  research. 

The  basic  concepts,  however,  may  be  extended  to  any  geometry 
and  any  flow  condition,  provided  of  course  that  these  differ- 
ences are  taken  into  account  in  interpreting  observed  dis- 
charge characteristics. 

A.  BREAKDOWN  VOLTAGE  (VL) 

b 

1 . Uniform  Field  (parallel  plate) 

Cobine  (Ref.  7)  and  von  Engle  (Ref.  8)  both  analyze 
the  effects  of  space  charge  on  the  current  and  voltage  dis- 
tributions in  the  gap  between  parallel  plate  electrodes. 

The  case  of  the  positive  ion  emitting  anode  from  Reference 
(8)  is  suitable  for  the  purposes  of  the  present  discussion. 
Consider  first  a uniform  electric  field  between  the  plates 
with  no  space  charge.  The  electric  field  will  be  constant 
from  anode  to  cathode  at: 

E=V/d 

where  V=applied  voltage 

and  d=electrode  spacing. 

Now,  consider  the  effect  of  a space  charge  modifica- 
tion of  the  field  structure,  with  the  anode  emitting  posi- 
tive ions.  Assuming  that  the  bulk  of  the  discharge  is 
electrically  neutral  as  is  normally  the  case,  there  will 
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be  a small  region  near  the  anode  in  which  the  ion  density 
exceeds  the  neutral  charge  density.  A net  positive  charge 
will  then  characterize  this  area. 

The  effect  of  the  positive  space  charge  is  to  lower 
the  electric  field  in  this  region,  resulting  in  the  type 
of  field  and  voltage  distributions  shown  in  Figure  6.  The 
modified  E and  V curves  assume  constant  charge  density  out 
to  X*Xgc,  which  is  not  physically  possible,  but  it  illus- 
trates the  point.  It  should  be  noted  from  the  figure  that 
an  intense  space  charge  can  actually  cause  a positive 
initial  voltage  gradient.  At  point  P in  the  figure,  the 
voltage  gradient  is  zero,  and  point  P is  referred  to  as  a 
"virtual"  anode. 

Referring  to  the  distance  at  which  the  virtual  anode 
lies  from  the  origin  as  5 , it  can  be  seen  that  the  value 
of  the  electric  field  in  the  neutral  region  must  increase 
to  account  for  the  same  voltage  drop  between  the  electrodes. 
One  way  of  interpreting  this  is  to  define  an  equivalent  gap 
length 

d'-d-fi  (3) 

Noting  also  that  the  effective  voltage  drop  from  the 
virtual  anode  to  the  cathode  is  increased  by  the  amount  AV, 
it  is  possible  to  write  an  expression  for  the  effective 
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electric  field  strength  in  the  neutral  region  as: 

E.  _ V+  AV 
d' 

The  important  point  here  is  that  the  apparent  physical  gap 
length  is  reduced  by  space  charge  modification  of  the 
electric  field  structure  (Ref.  9).  This  "virtual"  electrode 
effect  applies  to  any  electrode  near  which  charges  of  the 
same  polarity  accumulate  (Ref.  7). 

In  a normal  parallel  plate  gap,  the  situation  is 
very  similar  to  the  one  just  described,  since  positive 
charge  tends  to  build  up  around  the  anode  (Ref.  4).  If 
this  charge  reaches  a high  enough  concentration  just  prior 
to  breakdown,  the  result  will  be  a "virtual"  anode,  and  the 
effective  gap  length  d'  will  reflect  this  reduction. 

Neglecting  AV,  then,  the  breakdown  field  strength 

would  be 

± Vb/d'  = Vb/(d-6  ) 
which  can  be  written  as 

Eb  = 6/d> 

and  for  <5/d « this  can  be  written  as: 

Eb  i <Vd)(1  + 5/d)  (*) 

This  equation  provides  a basis  for  a test  of  the  space-charge 
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model  concept  which  will  be  used  to  describe  processes  in 
the  non-uniform  field,  so  it  is  worthwhile  to  see  if  it 
gives  results  approximately  in  agreement  with  observations. 

Since  it  was  assumed  that  5/d  was  small,  E^  can 
be  evaluated  at  as  large  a value  of  d as  possible  to  mini- 
mize the  contribution  of  the  <5/d  term.  From  the  data  of 
Reference  (4),  at  d=20  cm,  V^=510  kV,  so 

= Vb/d)*>  & 510/20  = 25.5  kV/cm 

Now,  to  estimate  <5,  the  data  point  at  d=10  cm  can  be  used, 
where  the  measured  V^=265  kV.  The  hypothesis  is  that  for 
any  gap, 

Eb  = Vb/d)~  (1  + 6/d  = 265/10 

so,  for  the  10  cm  gap, 

(1  + 5/10)  = 26.5/25.5 

and 

6 = .39  cm 

Now,  the  final  result  is: 

^ \ «o(1  + 6/d)  = 25'5  (1  + *39/d)  <5) 

In  reality,  the  anode  extension  can  never  be  larger  than 
the  physical  gap,  so  the  upper  limit  of  E^  based  on  this 
crude  approximation,  is  2 E,  ^ = 51kV. 
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It  is  also  noted  that  in  the  absence  of  space  charge 


- 


effects  (i.e.  6 =0),  E^=const  which  is  the  required  result. 

It  can  be  seen  in  Figure  (7)  that  at  2 cm  and  above, 
the  computed  values  of  the  breakdown  field  are  quite  close 
to  the  measured  values;  closer,  in  fact,  than  the  values 
computed  from  the  other  models  listed  in  Reference  (4) . 

This  agreement  demonstrates  the  possible  utility  of  such  a 
space-charge  model,  and  this  concept  will  be  used  in  a 
detailed  analysis  of  gasdynamic  effects  on  a non-uniform 
field. 

2 . Non-Uniform  Fields 

The  basic  concept  described  above  for  the  uniform 
field  can  also  be  applied  to  the  non-uniform  field  distri- 
bution used  in  the  present  experimental  investigation. 

First,  the  process  will  be  described  as  in  a non-flowing 
discharge  section,  and  later,  the  effects  of  gas  convection 
and  turbulence  will  be  discussed.  Experimental  data  will 
then  be  scrutinized  in  light  of  the  arguments  presented. 

a.  No -Flow 

First,  consider  a stationary  discharge  in  a 
point-plane  gap.  For  sharp  points,  the  region  near  the 
point  is  characterized  by  very  high  electric  fields  just 
prior  to  spark  breakdown  of  the  gap.  For  a simple  model, 
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consider  the  field  distribution  for  a hyperboloidal  point- 
to-plane  (Ref.  10): 


E(x) 


E 


P 


2V 

R log(4  d/R) 


(6) 


Ep  - field  at  surface  of  point 
R ■ point  radius 
V - applied  voltage 
d = gap  length 

For  typical  values  of  R and  d of  0.05  and  3 cm 
respectively,  breakdown  voltage  has  been  found  experimentally 
to  be  about  27,400  volts  in  atmospheric  air.  This  results 
in  an  Ep  of: 

E = 2(27,40,0^ ± 200,000  v/cm 

P 0.05  log(4*3/.05) 

This  corresponds  to  a reduced  field  (E^/P)  of  263  v/cm-torr. 

Since  this  value  of  electric  field  is  well  above 
the  atmospheric  air  ionization  field  strength  of  15,000 
v/cm  (20  v/cm-torr)  (Ref.  11),  it  is  evident  that  intense 
ionization  of  the  air  will  be  taking  place  in  the  vicinity 
of  the  point  by  high  energy  electrons  accelerating  toward 
the  positive  point. 
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As  an  electron  accelerates  toward  the  point, 
it  creates  an  electron  avalanche.  Since  the  electron 
mobility  is  much  higher  than  the  ion  mobility,  the  elec- 
trons speed  into  the  positive  point  and  get  absorbed 
leaving  behind  the  ion  space  charge  in  front  of  the  pin 
tip.  Detailed  descriptions  of  the  development  of  streamers 
and  the  electron  avalanche  can  be  found  in  References  (7, 

8,  10,  12).  It  is  sufficient  for  the  present  to  realize 
that  dense  positive  space  charge  can  be  formed  in  the 
immediate  vicinity  of  the  anode  which  can  significantly 
alter  the  electric  field  structure  there.  (In  an  electro- 
negative gas  such  as  air,  negative  ions  would  also  be 
formed,  and  their  number,  to  a first  approximation,  can  be 
regarded  as  a fraction  of  the  number  of  positive  ions.) 

The  net  charge  density  in  the  vicinity  of  the  anode,  there- 
fore, remains  positive,  and  negative  ion  effects  will  be 
disregarded  for  the  time  being  (Ref.  9).  The  positive  ions, 

under  the  influence  of  electrostatic  repulsion  and  diffusion, 

* 

begin  to  move  slowly  away  from  the  positive  point,  creating 
the  space  charge  region.  As  voltage  is  raised,  the  dis- 
charge passes  through  a stable,  pulseless  "glow"  mode,  and 
then,  as  applied  voltage  approaches  the  gap  breakdown  value, 
the  formation  of  pre-breakdown  streamers  causes  the  frequency 


of  current  pulses  to  increase  to  approximately  10  to 
10^  Hz.  The  maximum  length  of  streamers  measured  under 
pulsed  conditions  has  been  observed  to  be  a linear  function 
of  the  applied  voltage  when  no  space  charge  exists  to  dis- 
tort the  electric  field  (Ref.  4).  It  seems  reasonable, 
then,  that  the  pre-breakdown  streamers  formed  in  a static, 
steady-state  discharge  would  behave  in  a similar  manner. 

The  high  concentrations  of  positive  charge  in  the  vicinity 
of  the  anode  cause  an  apparent  anode  extension  into  the  gap, 
reducing  the  effective  gap  length  much  as  in  the  case  of  the 
"virtual"  anode  described  in  the  previous  section.  These 
pre-breakdown  streamers  have  only  to  travel  the  reduced, 
effective  gap,  then,  to  effect  breakdown  of  the  gap  by  the 
processes  described  in  Reference  (4) . If  this  assumption 
is  correct,  then  influencing  the  extent  of  the  space -charge 
region  by  flow  or  turbulence  should  show  up  as  changes  in 
the  measured  values  of  breakdown  voltage.  To  test  this 
assumption,  it  is  necessary  to  estimate  what  the  original 
size  of  the  space  charge  region  is  under  no-flow  conditions. 
The  following  sections  address  themselves  to  this  problem, 
and  then  to  the  problem  of  determining  the  changes  in 
breakdown  voltage  which  might  be  expected  due  to  alteration 
of  the  space-charge  region  by  gasdynamic  processes. 
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It  is  interesting,  and  will  be  useful  later,  to  examine 
the  effects  of  gas  density  variation  on  the  anode  exten- 
sion and  breakdown  voltage  under  no-flow  conditions.  The 
following  relationships  are  used: 


Mean  Electron  Energy 
for  Ionization: 

eE  A eV^ 

(7) 

Electric  Field 
Distribution: 

E(x)  • (R^)EP 

(8) 

Breakdown  Voltage 
Relation : 

Vb  * ad'+b 

(9) 

Virtual  Anode  Effect: 

d'  = d-  6 

(10) 

Mean  Free  Path: 

A = 4V2~  A = -■  - — - 

e 5 P(torr) 

(ID 

e = electronic  charge  = 1.6x10"  coul 

= ionization  potential  = 15.5  volts  for  air 
a,b,c  = constants 

A = electron  and  gas  molecule  mean  free  paths 
e > § 

P = pressure 

and  all  other  symbols  have  been  previously  defined. 

Solving  for  the  ionization  field  strength 

from  Eq . ( 7)  , 

Ei  * Vt/  Xe 


and  using  Eq.  (11),  this  becomes 


Ei  * ViP/c 


Now,  using  this  value  in  Eq.  (8),  it  is  possible  to  solve 
for  a characteristic  "ionization  distance"  (the  maximum 
distance  from  the  anode  at  which  a free  electron  would,  on 
the  average,  in  one  mean  free  path,  gain  sufficient  energy 
to  ionize  by  collision).  Rearranging  Eq.  (8)  with  E(x)  = E^, 
and  recalling  E^  from  Eq.  (6): 

(R+x.)  - — — E = R EP(v)  = -V. 
i Ei  p (V-P/c)  P 

It  is  assumed  that  the  anode  extension  6 is  proportional 
to  x^,  and  that  R<<6  , so  that 


5 i x = _v_  i v_  (Const  temp) 
* P p (Ideal  gas) 


Now,  in  direct  analogy  with  the  uniform  field  example  of 
the  previous  section,  the  effective  gap  length  is  expressed 


by  Eq.  (10)  as 


d'  = d-6  = d-  — c' 


and  it  is  seen  that  the  effective  gap  length  is  reduced  as 

the  gas  density  is  reduced.  If  this  value  of  d'  is  inserted 

into  Eq.  (9) , after  some  rearrangement  it  can  be  shown  that 

V = .Cgd+b)  p (13) 

b ( P+  ac') 
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The  form  of  this  equation  will  be  shown  to  describe  the 
actual  breakdown  voltage  density  dependence  fairly  well. 
For  example,  if  it  is  rewritten  as 


where 


VVbo 


2(  p / Pn) 

( P/  P0+l) 


(14) 


V^Q  = ad+b  = at  one  atmosphere 

the  values  of  V^/V^  vs.  P / P0  agree  with  the  observation. 

The  constant  term  in  Eq.  (14)  is  dependent  on  the  specific 
electrode  geometry  and  type  of  gas.  Table  II  is  a collec- 
tion of  data  for  a positive  point-to-plane  for  gap  lengths 
and  densities  in  the  present  range  of  interest. 

In  making  comparisons  of  measured  data  at 
one  atmosphere  pressure,  it  was  noted  that  Bandel’s  data 
from  Table  II  for  the  8 cm  gap  is  inconsistent  with  the 
data  of  References  (6)  and  (7)  for  a point-plane  gap  of 
8 cm,  and  with  trends  noted  in  the  present  experimental  work. 
Also,  due  to  the  fact  that  Bandel's  breakdown  voltage  at  one 
atmosphere  was  nearly  the  same  for  4 and  8 cm  gaps,  there 
is  doubt  that  the  value  of  given  is  appropriate  for  the 
present  analysis.  Therefore,  the  value  of  used  in  the 
calculation  of  the  voltage  ratio  for  the  8 cm  gap  for  com- 
parison with  Eq . (14)  was  60  kV,  which  is  in  close  agreement 
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with  the  other  references  cited.  Bandel's  second  point 
for  the  8 cm  gap  was  also  discarded  for  the  reasons  already 
mentioned.  Figure  (8)  shows  the  results  of  observed  values 
of  the  voltage  ratio  vs.  the  predicted  values  of  Eq.  (14). 
The  close  agreement  tends  to  support  the  "virtual"  anode/ 
space-charge  effects  model  outlined  above. 


TABLE  II 


Density 

Ratio 

Vb 

(kV) 

*b  < 

a) 

Meas 

/ Calc 

Meas  / 

Calc 

Gap 

Length 

Source 

.05 

2.0 

2.6 

0 

0 

4 

Bandel 

.09 

4.0 

5.8 

65 

60 

4 

(Ref  18) 

ff 

.15 

7.2 

8.7 

105 

110 

4 

ft 

.24 

14.2 

14.7 

225 

180 

4 

It 

.39 

22 

21 

240 

230 

4 

1! 

.63 

33 

33 

240 

230 

4 

ft 

.98 

38 

38 

130 

130 

4 

If 

.05 

4.0 

5.4 

25 

20 

8 

ff 

.09 

8.0 

9.0 

150 

90 

8 

ft 

.15 

15 

14.3 

180 

166 

8 

ft 

.24 

24 

21.8 

215 

204 

8 

ft 

.39 

32 

31 

170 

185 

8 

fl 
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TABLE  II  (Continued) 


Density 

Ratio 

Vb 

(kV) 

Xb 

( a) 

Meas 

/ Calc 

Meas 

/ Calc 

Gap 

Length 

Source 

.63 

(33) 

40 

105 

105 

8 

Bandel 

.98 

(38) 

41 

40 

36 

8 

(Ref  18) 

II 

1.0 

4.8 

2.2 

8.5 

1.7 

.3 

Werner 

1.0 

5.0 

5.8 

5.5 

4.8 

.6 

(Ref  15) 

II 

1.0 

9.2 

11 

14 

10.6 

1.0 

Townsend 

1.0 

10 

18 

(10) 

21.5 

1.5 

(Ref  17) 

II 

1.0 

25 

25 

• mm 

37.5 

2.0 

Meek  & Cragg; 

1.0 

60 

-- 

-- 

-- 

8.0 

(Ref  10) 

II 

1.0 

65 

• — 

- — 

- — 

8.0 

von  Engle 

(Ref  7) 

(See  Appendix  B for  calculation  procedures  and  formulae.) 

(NOTE:  Data  in  parenthesis  do  not  seem  to  correlate  with 

trends  observed  in  other  references.) 

b.  Gas  Convection/Flow  Orientation  Effects  Analysis 
It  has  been  observed  from  our  work  that  the 
orientation  of  the  flow  field  relative  to  the  electric  field 
has  a marked  effect  on  the  electrical  properties  of  the 
discharge.  Most  evident  is  the  change  of  breakdown  voltage 
observed  in  the  various  orientations  tested  (Fig.  9). 
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Figures  (10,  11)  illustrate  the  nature  of  these  changes. 

It  can  be  seen  from  the  Figures  that  convection  increases 
Vfa/Vbo  in  the  reverse  and  cross-flow  modes  (Fig.  11) 
quickly  to  a nearly  constant  value  on  the  order  of  1.25. 

In  the  ordinary  mode,  however,  the  general  trend  is  for  an 
initial  decrease  in  v^/v^0  to  a minimum  in  the  neighborhood 
of  0.9-0.95  at  around  75  m/sec  flow  velocity.  After  that, 
the  ratio  begins  to  increase  again,  occasionally  to  a 
value  greater  than  one. 

A space-charge-convection-effects  model  has 
been  developed  which  successfully  explains  some  of  the 
major  features  of  the  breakdown  voltage  behavior. 

Two  orientations  were  tested,  one  longitudinal 
mode  with  flow  from  anode  to  cathode  (ordinary  mode)  (Fig. 
9a),  and  one  cross-flow  mode  with  the  discharge  being 
transverse  to  the  gas  flow  (Fig.  9c).  The  third  mode 
analyzed  by  the  present  model  is  another  longitudinal  mode, 
as  in  case  one,  but  with  the  flow  direction  reversed  so 
that  the  gas  flow  is  from  cathode  to  anode  (Fig.  9b) . 

(1)  Space-Charge-Convection-Effects  Model . 
Although  the  discharge  sections  used  consisted  of  from  33- 
39  positive  pins  spaced  one  centimeter  apart  for  the  anode, 
and  a grounded  plane- type  cathode,  the  electric  field 
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structure  for  any  given  pin  is  nearly  identical  with  that 
of  a single  point-plane  configuration,  since  the  pins  are 
far  enough  apart  that  the  space  charge  regions  do  not  over- 
lap (Fig. *2).  Therefore,  the  model  is  discussed  in  terms 
of  the  dynamic  processes  taking  place  at  one  pin,  with  the 
understanding  that  the  processes  apply  equally  to  the  test 
section  as  a whole,  provided  that  all  pins  contribute 
uniformly  to  the  current. 

It  must  be  noted  at  this  point,  however, 
that  due  to  the  statistical  nature  of  breakdown  streamer 
size  and  intensity,  breakdown  becomes  a probability  related 
process.  For  example,  if  the  voltage  across  a single  point 
gap  is  set  so  as  to  give  a probability  of  breakdown  (Pfp 
in  a certain  time  interval  (At)  of  Pf^=.l,  then  its  prob- 
ability of  not  breaking  down  is 


psi=1-pfr-9 


(15) 


Now,  if  there  are  10  pins  acting  inde- 
pendently at  the  same  applied  voltage,  each  with 
then  the  probability  of  the  10  pin  gap  not  breaking  down 
would  be  given  by  the  following  expression: 


PsKf  l-(lp£i>10  = l-(l-.l)10  = .65  (16) 

This  demonstrates  the  increased  susceptibility  of  a multi-pin 
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configuration  (with  statistically  independent  pins)  to 
break  down.  Therefore,  to  obtain  the  same  breakdown  prob- 
ability, a multi-pin  gap  must  be  at  a lower  voltage,  and 
observed  values  of  breakdown  voltage  in  multiple  pin  con- 
figurations may  reflect  this  reduction. 

The  expression  used  for  the  electric  field 
in  the  space-charge  convection  model  is  the  same  as  used 
in  previous  sections  i.e.,  that  for  the  hyperboloidal  point- 
to-plane.  In  this  model  the  positive  space  charge  is 
assumed  to  be  at  a high  enough  density  to  cause  significant 
distortion  of  the  undisturbed  field  out  to  a distance  where 
the  undistorted  positive  point  field  strength  is  = 20,000 
volts/cm.  This  point  in  the  field  was  chosen  because  it 
is  well  inside  the  region  of  ionization  by  electron  collision 
of  = 15,000  volts/cm  for  air  at  atmospheric  pressure.  The 
corresponding  reduced  field  strength  (E/P)  - 20  v/cm-torr. 
Referring  again  to  Figure  (1) , an  electron  arriving  at  point 
a,  E = 15,000  v/cm,  would,  on  the  average,  have  enough 
energy  to  ionize  air  molecules  by  collision  and  electron 
avalanching  would  start  toward  the  positive  tip.  The 
result,  as  described  in  the  previous  section,  is  positive 
space  charge  accumulation  in  front  of  the  pin.  The  posi- 
tion where  the  pin  field  strength  approaches  20,000  v/cm, 
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point  b,  is  conjectured  to  be  in  the  critical  region  where 
the  space  charge  field  becomes  significant  compared  with 
the  undistorted  field.  When  this  happens,  the  rest  of  the 
gap  acts  as  if  the  anode  lies  inside  the  high  field  region 


caused  by  the  dense  positive  space  charge,  and  the  gap  is 
effectively  shortened.  Figure  (12)  is  an  illustration  of 
this  idea. 

The  distance  defined  above  is  used  in  the 


I 

rest  of  the  analysis  as  the  magnitude  of  the  anode  extension, 
previously  defined  as  5 in  the  no -flow  analysis  of  the 
density  dependence  of  the  breakdown  voltage.  The  effective 
gap  length,  as  proposed  by  the  model,  would  be  the  distance 
between  the  cathode  and  the  point  at  which  the  static  elec- 
tric field  first  reaches  20,000  v/cm.  Expressed  mathematically, 
this  would  be 


d'  = d - 


620K 


where  the  subscript  indicates  the  field  value  at  which 
is  computed.  This  value  of  6 will  represent  the  no-flow 
anode  extension  and  will  be  referred  to  as  <50=  SJ0K  . 

(2)  Experimental  Inputs  and  Basic  Model. 
Experimental  data  over  a range  of  gap  lengths  under  no-flow 
conditions  at  atmospheric  pressure  have  yielded  an  approxi- 
mately linear  relation  between  gap  length  and  breakdown 
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voltage,  closely  approximating  the  point-to-plane  curves 
of  References  (8,  10)  (Fig.  13). 

The  apparent  reduction  in  gap  length, 
can  again  be  calculated  from  Eq.  (8),  using  E(x)  = 20,000 
v/cm,  as 


6 


o 


X20K 


REp 

20,000 


-R 


(17) 


but  note  that 

Ep  = 2V 

R log(4d/R) 

Using  the  experimentally  derived  relation  for  breakdown 
voltage  (Fig.  13) : 

Vb  = 7800d+4000  (volts)  (18) 

With  an  estimated  R=0.05  cm,  then 


<5  (d) 


. 05(  2(  7800d  + 4000)) 
.051og(4d/ 105) 
20,000 


(19) 

-.05 


TABLE  III 


d(cm) 

<5  (cm) 

a /d 

2 

.34 

.17 

4 

.56 

.14 

The  quantity  6/d  in  Table  3 is  a measure 
of  the  relative  reduction  in  gap  length  for  two  gaps  in 
the  range  of  interest.  It  can  be  seen  that  the  relative 
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reduction  due  to  the  presence  of  space  charge  over  the 
2 cm  difference  is  only  about  37<> . What  this  means  is  that 
the  linearity  suggested  in  the  data  may  not  be  between 
and  d,  the  actual  physical  gap  length,  but  may  be  a linearity 
between  Vb  and  d' , the  apparent  gap  as  speculated  earlier. 

Due  to  the  small  relative  reduction  of  gap  length,  it  is 
not  possible  to  determine  from  the  data  alone  which  is  the 
actual  case. 

In  the  ensuing  discussion,  it  is  assumed 
that  is  a linear  function  of  the  effective,  reduced  gap: 

Vb  i d'  = d - 6 (9,  10) 

Adoption  of  this  hypothesis  leads  to  some  interesting 
conclusions . 

( 3)  Calculation  of  Maximum  Voltage  Ratio  due 
to  Convection.  In  flow  orientations  two  and  three,  (Fig. 

9 b&c) , it  can  be  deduced  from  the  above  theory,  that  the 
maximum  value  of  breakdown  voltage  occurs  when  convection 
effects  remove  all  the  space  charge  from  the  gap,  thereby 
increasing  the  effective  gap  spacing  from  its  normal  no- 
flow value  of  d^  = (d  - 6q)  to  d,  its  upper  limit.  The 
maximum  ratio  attainable  then  at  any  velocity  should  be 

<Vvbo>ma*  - -d  t~~  (20) 

O 

where  the  subscript  'o'  denotes  the  no-flow  value  (6„  = <5,_,.  ). 

v JLO  ^ 


The  Table  below  summarizes  the  results  of 


this  calculation  for  several  gaps  in  the  range  of  the  ex- 
perimental data. 

TABLE  IV 

REVERSE  AND  CROSS -FLOW  MAXIMUM  VOLTAGE  RATIOS 


Gap  Length 
d(  cm) 

Anode  Ext 
(Eq.  17) 
6o  (cm) 

Effective  Gap  Voltage  Ratio  (max) 

do=d-  50(cm)  Theory  Experiment 

1 

.22 

.78 

1.28 

-- 

1.4 

.27 

1.13 

1.24 

1.27 

1.45 

.27 

1.18 

1.23 

1.26 

2 

.34 

1.66 

1.20 

— 

3 

.45 

2.55 

1.18 

1.20 

3.2 

.47 

2.73 

1.17 

1.15 

4 

.56 

3.44 

1.16 

1.13 

(See  also  Figure  14) 

It  is  evident  from* the  Table  that  theory 
and  experiment  agree  quite  well,  despite  many  crude  approxi- 
mations which  went  into  the  calculations.  The  close  numeri- 
cal agreement  tends  to  verify  the  initial  assumption  that 
the  point  where  E(x)=20,000  v/cm  may  be  considered  to  be 
the  approximate  axial  range  of  the  anode  extension  under 


no-flow  conditions.  Perhaps  even  more  significant  than 


the  numerical  agreement  is  the  fact  that  the  theory  pre- 
dicts a decrease  in  the  maximum  voltage  ratio  with  increas- 
ing gap  length.  This  trend  has  been  experimentally  verified. 


(4)  Minimum  Voltage  Ratio.  The  ordinary  flow 
mode,  Figure  (9a),  is  more  complicated  as  can  be  seen  from 
the  data  summary  plot  of  Figure  (10). 

It  might  be  expected  from  a space  charge 
convection  point  of  view,  that  in  this  orientation  the 
space  charge  would  be  carried  toward  the  cathode  in  some 
characteristic  time  interval  which  would  tend  to  increase 
<5  . If  6 increases,  the  effective  gap  length  decreases, 
which  in  turn  results  in  a decrease  in  V^,  since  - d ' . 
This  agrees  with  the  initial  trend  of  the  data  of  Figure 
(10).  The  breakdown  voltage  vs.  velocity  curves  show  a 
minimum  though,  beyond  which  begins  to  increase  again. 

At  short  gaps,  in  fact,  has  been  observed  to  not  decrease 
at  all,  but  show  a slight  increase  which  remains  relatively 
constant  with  increasing  velocity. 

In  order  to  shed  some  light  on  these 
apparent  inconsistencies  at  short  gaps,  and  to  attempt  to 
explain  the  complex  nature  of  the  breakdown  in  this  orienta- 
tion, it  is  necessary  to  look  into  the  relative  magnitudes 
of  the  convective  and  field  effects  on  the  ions  with 
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reference  to  their  specific  locations  in  the  total  field 
structure.  A convenient  parameter  for  such  a comparison 
is  the  "slip  parameter" 


g _ _U_  = U_  * gas  velocity  _ convection  effect  (21) 
kE  v^  drift  velocity  electrostatic  effect 

When  the  order  of  magnitude  of  the  slip 
o 

parameter  is  = 1,  convection  effects  become  as  important 
to  ion  motion  and  location  as  the  E field,  whereas  if 
S = U/v^<<  1,  convection  effects  would  be  negligible. 
Conversely,  for  S»  1,  the  electrostatic  effects  become 
negligible.  Figures  (15  & 16)  illustrate  the  magnitude 
of  this  parameter  plotted  against  U for  the  field  point 
E=20,000  v/cm,  and  E^Ep. 

It  is  evident  from  Figure  (15)  that  an 
ion  located  in  an  electric  field  of  strength  E=20,000  v/cm, 
would  experience  significant  convective  effects  at  very 
low  velocities.  Figure  (16)  is  an  illustration  of  the 
slip  parameter  values  at  the  anode  surface,  S=U/v^,  with- 
out considering  the  effects  of  space  charge.  However,  it 
must  be  remembered  that  in  actuality,  a positive  ion  at 
the  anode  surface  is  located  between  the  tip  and  a region 
of  high  positive  space  charge,  which  reduces  the  effective 
field  strength  at  its  location  (See  Ep'  in  Fig.  12). 
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By  reducing  the  effective  value  of  Ep  by 
an  amount  f(n"*")  , then  E ' =E  -f (n+) , and  S =U/k+E  ' increases, 
which  could  render  convection  a significant  effect  even  at 
the  tips  of  the  pins.  An  important  aspect  of  the  quantity 
U/k+Ep'  is  that  it  is  dependent  on  the  gap  length,  since, 
as  stated  earlier,  at  breakdown: 

= 7300d  + 4000  (volts)  and 

E\  = - f(n+)  which  yields 

Pb  Rlog(4d/R) 

U U 

vdp'  k't‘2(7800d+4000)  /Rlog(4d/R)  -f  (n'r) 

Due  to  the  nature  of  the  log  term,  this  yields  approximately 
the  following  important  dependence: 


U 

d-fCn*-) 


jU/dl=_U  (l+f>iU  f i+c(i)l 
1-f  d d d [ d 


(22) 


What  this  equation  effectively  says  is  that 
at  short  gaps,  the  convective  effects  will  become  more 
prominant  and  may  tend  to  dominate  the  discharge,  while  at 
longer  gaps,  at  least  through  the  range  tested,  both  factors 
may  be  important.  It  is  noteworthy  that  the  quantity  f(n+) 
is  proportional  to  the  current,  since  the  charged  particle 
density  is  proportional  to  the  current.  This  factor  would 
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also  tend  to  increase  the  slip  parameter  as  the  current 
increased. 

The  alternative  to  attempting  to  analyse 
the  exact  magnitude  of  the  effective  electric  field  at  the 
anode  surface  is  to  inspect  the  available  data  for  evidence 
of  significant  flow  effects  which  may  alter  the  breakdown 
voltage  characteristics  due  to  creation  of  a sufficiently 
high  slip  parameter  at  the  anode  surface. 

A closer  inspection  of  the  data  of  Figure 
(10)  indicates  a reversal  of  the  declining  trend  in  vs.  U 
for  those  gaps  greater  than  about  2 cm  in  the  range  of  60- 
90  m/sec  flow  velocity.  Therefore,  it  is  assumed  that 
this  is  the  velocity  region  in  which  convection  begins  to 
seriously  affect  the  ion  population  closest  to  the  pins. 

The  current  vs.  velocity  curve  also  indicates  that  the 
current  starts  to  rise  sharply  in  this  range,  so  the  com- 
bined effects  of  high  velocity,  high  current,  and  short 
gap  (Eq.  22),  appear  to  significantly  alter  the  dynamics  of 
the  situation.  Since  convection  becomes  important,  it  is 
reasoned  that,  in  a certain  characteristic  time  interval 
(The  time  between  current  pulses  or  breakdown  streamers 
seems  to  be  appropriate  in  this  case,  and  can  be  estimated 
from  the  pulse  frequency  cited  in  the  introduction  as 
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an  amount  away  from  the  pin  tip  anodes,  that  the  space 
charge  no  longer  acts  as  an  anode  extension.  Instead  it 
becomes  a detached  region  of  floating  charge,  no  longer 
able  to  act  as  a bridge  for  infeeding  electrons  at  the 
head  of  succeeding  avalanches,  to  run  all  the  way  into  the 
anode  itself.  Figure  (17)  illustrates  the  various  stages 
of  progression  of  the  events  just  described.  A separate 
shortening  effect  of  the  anode  extension,  due  to  high 
current  density  and  recombination,  may  also  play  a part 
in  the  observed  voltage  characteristic  rise,  and  it  will 
be  discussed  fully  in  the  next  section  in  connection  with 
observations  in  the  cross-flow  mode. 

To  summarize  up  to  this  point  for  the 
ordinary  flow  mode,  the  sequence  of  events  as  described  is 
as  follows: 

(a)  No  Flow 

S = U/v^  = 0 Fig.  (17a) 

The  physical  gap  is  effectively  reduced 

by  the  concentration  of  space  charge.  The  effective  gap 

is  d'  = d-8  , and  V.  /V.  -1.0. 
o o b bo 

(b)  Intermediate  Flow  (IF) : 


S ' U/vd(20K)2  1 Fig‘  (17b) 
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Convective  effects  extend  8 further  into  . 


the  gap  and  = d^  -UT  . ^b^bo< 

(c)  Fast  Flow  (FF) : 

S - U/v  ' 2 1 Fig.  (17c) 

P dp 

Space  charge  detaches  from  tip  and 
d'pF > d ' i so  V^/V^o  begins  to  increase  again. 

Biblarz  and  Nelson  (Ref.  16)  stated  that 
the  maximum  observed  voltage  ratio  in  the  ordinary  flow 
mode  equaled  but  did  not  exceed  the  maximum  reverse  flow 
value  of  1.13  for  the  3.9  cm  gap  in  their  tests.  This 
observation  is  consistent  with  the  present  model,  since 
1.13  is  the  maximum  ratio  for  any  flow  mode  and  represents 
the  entire  physical  gap. 

There  still  remains  the  strange  behavior 
of  the  breakdown  voltage  curve  for  the  short  gap  depicted 
in  Figure  (10).  Going  back  to  Eq.  (22),  recall  that 
U/vdp  = -^-(1+c^),  so  it  would  be  expected  that  at  suffi- 
ciently short  gaps  and  high  currents,  convection  effects 
would  dominate  the  space  charge  at  the  tips,  even  at  very 
low  velocities.  This  would  lead  to  early  (low  velocity) 
separation  of  the  space  charge  from  the  anode,  resulting 
in  an  early  rise  in  breakdown  voltage  probably  to  a near 
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constant  value  at  or  below  (V,  /V,  ) , N since,  as 

b bo  (max) 

decreased,  d'— ►d.  The  data  for  the  1.6  cm  gap  appear  to 
be  a manifestation  of  this  phenomenon. 

A numerical  example  may  help  to  clarify 
the  ideas  outlined  above. 

Let  d=3  cm  gap  length 


Vb-7800d+4000  = 27400  volts  (18) 

Location  of  field  point  E=20,000  v/cm 
from  Eq . (17) : 

820K  * cm 

Effective  gap  length  (due  to  space  charge 
presence) : 

d ' = d - 80  = 3 - .45  = 2.55  cm 

Effective  electric  field  (neglecting  AV): 

E'=Vb/d^  = 27000/2.55  = 10,588  v/cm 

Assume  Tf = 2(10)  ^ sec,  U=  76  m/sec  = 

7600  cm/sec 

Anode  extension  due  to  convection: 
x = UTf  - 7600(2xl0~5)  = . 15  cm 
d^  = dg-uTf  = 2. 55-. 15  = 2.40  cm 
and  finally, 

VVbo=E'^/Vbo=10-588(2*4)/27-°  = — 

Experimental  values  for  Vb/Vbo  in  the  76  m/sec  velocity 
range  fall  in  the  .94-. 96  range,  with  the  current  relatively 
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constant  for  all  measurements.  Above  this  velocity,  con- 
vection at  the  tip,  and/or  the  decreasing  recombination 
distance  (to  be  explained  in  next  Section)  takes  over,  and 


1 


the  effective  gap  length  and  Vb  begin  to  increase  again. 

(5)  Turbulence  Effects  on  Breakdown  Voltage. 
Figure  (18)  is  a schematic  of  visually  observed  spatial 
characteristics  of  the  tip  glow  patterns  around  the  pins 
under  no-flow,  laminar  flow,  and  turbulent  flow  conditions. 
It  can  be  seen  that  the  intense  glow  region  of  (a)  is 
greatly  extended  by  convection,  (b) , and  that  the  addition 
of  turbulence  tends  to  spread  the  glow  over  a much  larger 
region.  In  fact,  turbulence  appears  to  shorten  the  maxi- 
mum extent  of  the  glow  under  some  conditions.  It  would 
seem  probable  that  this  spreading  and  reduction  of  the 
glow  region  would  have  an  effect  on  the  breakdown  voltage. 
In  fact,  according  to  the  model,  it  should  increase  V^, 
due  to  physical  shortening  of  6 . This  trend  is  well 
illustrated  in  Figure  (10),  where  the  curves  for  the  2.9  cm 
gap  show  a definite  upward  shift  with  increasing  turbulence 
level.  The  symbols,  in  order  of  increasing  turbulence, 
are  v , ° , . It  should  be  noted  that  these  symbols 

also  indicate  the  order  of  increasing  current,  so  again, 
increased  current  density  could  cause  faster  recombination 
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and  smaller  anode  extension  which  would  also  cause  an 
increase  in  V^.  Figures  (19,  20)  are  photographs  of  the 
discharge  section  which  illustrate  the  patterns  of  Figure 


(18). 


The  next  section  analyzes  in  more  detail 


the  nature  of  flow  and  charged  particle  interactions  in 
the  cross-flow  mode. 

(6)  Cross-Flow  Stability  and  Field  Interactions. 
Figure  (21)  is  a schematic  representation  of  the  space  charge, 
electric  and  cross-flow  field  interactions  as  envisioned 
with  the  present  model.  The  visual  character  of  the  dis- 


charge, as  observed  in  many  tests,  appears  to  bear  out  the 
conclusions  reached  in  the  following  analysis. 

Under  no-flow  conditions,  the  discharge 
gap  experiences  an  apparent  anode  extension  80,  due  to 
space  charge  accumulation  as  in  Figure  (17a).  As  pointed 
out  in  the  last  Section,  this  apparent  extension  results  in 
a reduction  of  effective  gap  length  to  d^=d-  6 . It  was 
assumed  that  the  breakdown  voltage  was  proportional  to  d^, 
so  the  maximum  voltage  ratio  attainable  in  this  configura- 
tion would  be  realized  if  the  space  charge  could  be  removed 
completely,  resulting  in  d^-*-d  and  V^/V^^-^d/(d-  6q)  . 
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Due  to  the  nearly  orthogonal  orientation 
of  the  electric  and  flow  fields  in  the  cross-flow  arrange- 
ment, these  two  effects  operate  nearly  independently  of 
each  other  to  impart  motion  to  any  ionic  space  charge  under 
their  respective  influences.  For  example,  a positive  ion 
formed  at  the  tip  of  pin  (a)  in  Figure  (21)  would  experience 
electrostatic  drift  and  convective  motion  as  follows: 

X direction:  vx*U 

Y direction:  Vy=k+E(y) 

The  important  result  is  that  the  space  charge  location  is 
strongly  influenced  by  the  magnitude  of  U.  Taking  the  tip 
of  pin  (a)  as  the  origin,  the  downstream  displacement  is: 

s = Ut  (cm) 

U = gas  velocity  (cm/sec) 
t * Tf  = time  (sec) 

For  U = 30  m /sec  (3000  cm/sec),  and 
t =*  Tf  ■ 2x10  ^ sec,  the  displacement  x = 0.06  cm.  For 
small  tip  radii,  this  amount  of  displacement  can  severely 
alter  the  local  field  structure,  resulting  in  an  asymetric 
field  distribution.  Consider  the  region  between  the  dis- 
placed space  charge  and  the  anode  for  example;  this  region 
would  be  characterized  by  a relatively  weak  field  as  the 
opposing  positive  point  and  space  charge  fields  tend  to 
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cancel  each  other.  On  the  other  hand,  the  region  from 
which  the  space  charge  was  removed  will  be  at  near  the 
full  undisturbed  point  potential.  The  undistorted  high 
field  region  on  the  upstream  side  of  the  pins  would  then 
be  a source  of  continuous,  stable  corona,  as  electron 
avalanches  feed  into  it  and  the  space  charge  is  quickly 
removed  by  convection.  Visual  observations  of  a cross- 
flow  discharge  give  the  impression  that  the  stable,  blue- 
violet  glow,  which  characterizes  region  1 of  Figure  (21), 
leans  into  the  wind.  This  result  at  first  seems  contrary 
to  intuition,  but  on  closer  inspection,  it  appears  to  be 
a natural  result  of  the  field  interactions. 

On  the  other  hand,  while  removal  of  space 
charge  from  region  1 enhances  stability  there,  convection 
of  that  charge  into  region  2 promotes  instability. 

It  has  been  shown  elsewhere  (Ref.  7)  that 
the  presence  of  space  charge  enhances  ionization  in  a uni- 
form gap  if  the  reduced  field  (E/P)  is  less  than  one  half 

* 

the  Stoletow  constant  B,  which  for  air  is  365  volts/cm-torr , 
and  B/2=183  v/cm-torr.  This  is  much  higher  than  E/P  in  the 
non-uniform  gap  used  in  the  present  investigation,  except 
in  the  immediate  vicinity  of  the  anode.  It  is  evident  then, 
that  convection  of  space  charge  into  region  2 will  result 
in  ionization  instabilities  there. 
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Manifestation  of  the  instability  described 
above  takes  the  form  of  numerous  arcs  across  the  gap,  com- 
mencing generally  from  the  second  and  third  rows  of  pins, 
and  spreading  to  the  first  row  of  pins  just  prior  to  break- 
down. Even  while  heavy  arcing  is  taking  place  from  the 
first  row  of  pins  however,  the  arcs  are  confined  to  the 
region  behind  the  blue-violet  glow  which  marks  the  beginning 
of  the  discharge  activity,  thus  demonstrating  the  destabi- 
lizing effect  of  the  space  charge. 

Since  the  power  supply  used  in  these  tests 
is  designed  to  "trip"  on  overvoltage  and  overcurrent  (and 
it  does  faithfully  in  the  ordinary  flow  mode) , the  numerous 
arcs  observed  must  be  limited  to  fairly  low  currents,  and 
it  is  not  until  applied  voltage  is  raised  to  the  value 
required  to  break  down  the  relatively  space-charge-free 
gap  across  the  first  row  of  pins  that  complete  breakdown 
occurs,  resulting  in  a large  current  pulse  with  resultant 
power  supply  shutdown. 

Actual  breakdown  streamer  propagation 

7 8 

velocities  have  been  reportedly  measured  in  the  10 7 to  10 
cm/sec  range  (Refs.  4,  10,  12),  so  that  the  breakdown 
streamer  can  travel  the  width  of  the  gap  (in  the  ranges 
tested)  virtually  unaffected  by  the  gas  velocity.  For 


example,  in  the  longest  cross-flow  g<*p  tested,  d=3.2  cm, 
the  time  of  transit  of  a primary  breakdown  streamer  is 
t=d/v  = 3.2/10^  = 3.2x10"^  sec.  The  displacement  of  an 
average  molecule  or  ion  by  convection  in  this  time  inter- 
val,  even  at  U=  150  m/sec  is  Ut  = 5x10  cm.  The  result 
is  that  in  the  cross-flow  mode,  alteration  of  the  field 
structure  begins  at  low  gas  velocities  by  space  charge 
convection,  resulting  in  a stable  upstream  region  of  very 
low  space  charge,  and  a downstream  region  of  instability 
caused  by  space  charge  distortion  effects.  Breakdown 
occurs  across  the  nearly  space-charge-free  gap  upstream 
at  a value  of  V^/V^  = d/(d-  5q)  , virtually  unaffected  by 
gas  convection  due  to  the  high  propagation  velocity  of  the 
breakdown  streamers. 

(7)  Observed  Effects  of  Charged  Particle 
Density  on  Recombination  Times.  A very  interesting  phenom- 
enon was  observed  during  the  course  of  the  cross-flow  tests. 
While  heavy  arcing  took  place  across  all  rows  of  pins,  un- 
distorted by  convection  due  to  their  high  propagation 
velocity,  a secondary  pattern  of  luminous  filamentary 
channels  was  also  evident  which  appeared  to  be  very  much 
affected  by  the  convective  velocity  (Fig.  21),  Contrary 
to  intuition  again,  however,  the  luminous  secondary 
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streamers  are  convected  a shorter  distance  at  higher 
velocities.  Their  intensity  also  increased  with  increasing 
velocity.  Figure  (22)  is  a graph  of  recombination  distance 
observations  vs.  velocity. 

In  order  to  understand  this  behavior,  it 
is  necessary  to  focus  on  the  exact  mechanism  which  is 
responsible  for  these  luminous  streamers.  In  the  streamer 
theory,  the  formation  of  a highly  ionized  channel  from 
anode  to  cathode  results  in  a highly  conductive  channel 
through  which  massive  quantities  of  secondary  electrons 
flow,  having  been  released  when  the  propagating  streamer 
strikes  the  cathode.  When  the  arc  is  extinguished,  a 
column  of  space  charge  remains  which  is  made  up  of  positive 
and  negative  ions  and  electrons.  When  radiative  recombina- 
tion takes  place  between  particles  of  opposite  charge, 
whatever  the  particular  specific  mechanisms,  energy  is 
released  in  the  form  of  photons.  The  secondary  streamers 
are  believed  to  be  caused  by  the  massive  recombination  of 
the  charged  particles  left  in  the  aftermath  of  an  arc. 

The  standard  approach  to  the  recombination 
problem  (Refs.  4,  7,  8",  13)  is  as  follows:  where  charges 

of  opposite  sign  reside  together,  they  will  be  mutually 
attracted  to  each  other,  and  so  will  come  together,  exchange 
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charge  and  neutralize  each  other,  emitting  radiation  in 
the  process.  There  are  several  varieties  of  recombination, 
all  of  which  are  discussed  in  detail  in  the  references 
cited.  It  is  sufficient  for  the  present  to  say  that  the 
recombination  rate  will  be  at  least  proportional  to  the 
numbers  of  charged  particles  available,  although  the  pro- 
portionality varies  with  the  process,  so  the  change  in  the 
number  of  existing  particles  in  a time  (dt)  will  be  given 
by  an  expression  of  the  form: 

dn  = - a n+n  dt  (23) 

r 

where  a ^ is  a constant  of  proportionality.  The  number  of 
negative  charges  is  not  always  equal  to  the  positive  charges, 
but  n can  be  expressed  as  a fraction  of  n+. 

n"  = cn+ 

where  c could  also  vary.  Now,  calling  n+  = n,  with  the 
aid  of  the  above,  Eq.  (23)  can  be  written  as 

dn  = - a cn^dt 
r 


and 


dn/n^  = - a cdt  (24) 

r 

which,  when  integrated  yields  the  following: 

1/n  - l/nQ  = <*rct 
showing  that: 


t 


rec 


= 1/  arc(l/n  - l/nQ) 


(25) 
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Reference  (13)  states  that  recombination  phenomena  are 
usually  recognized  by  the  linear  dependence  of  t on  1/n. 
But  n is  directly  proportional  to  j , the  current  density: 


j = e(n+v+  + n“v~)  = en(v+  + v"c) 


(26) 


and  also,  it  has  been  shown  by  other  experimenters  (Ref. 

14)  and  empirically  from  the  present  experiments  that 
2 

j & U - n.  Making  use  of  this  relation  between  n and  U 
in  Eq.  (25)  leads  to  the  following: 
tree  4 i/U2 

Multiplying  both  sides  by  the  gas  velocity  U gives  the 
distance  traveled  during  recombination  in  the  arc-caused 
space  charge  column: 

xrec  = ^Crec  = 

which  is  the  observed  dependence. 

It  is  possible  to  do  an  order  of  magnitude 
analysis  of  xrec  by  going  back  to  Eq.  (25).  Assuming  that 
c ° 1,  which  is  reasonable  under  these  circumstances,  and 

— O 

ar  = 2(10)  (cm  /sec)  (Ref.  13).  If  it  is  assumed  that 
the  final  number  of  charged  particles  is  = . lnQ  (907.  recom- 
bination) , then 

trec  = d/2*10'6)(l/.lno-l/no)  = (1/2*10_6)  (9/nQ)  (27) 

The  relation  between  n and  U can  be 

o 

estimated  by  solving  for  the  current  in  terms  of  n and  U 
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as  follows: 

j * en(k+E  + k_E)  (28) 

which  reduces  to 

j = en(k  E)  = env^  (29) 

since  k+«  k . The  drift  velocity  of  electrons  is  difficult 
to  estimate  due  to  the  variable  field  and  severe  space 
charge  modification,  but  a reasonable  approximation  is  the 
streamer  propagation  velocity  - 102  cm/sec  (Ref.  12) . 

This  gives 

j ° 1.6(10'19)(107)no=1.6(10'12)no  (30) 

Now,  from  observed  data,  j(U)  can  be 
calculated,  assuming  j ^ U . Taking  a point  from  the  data 
for  atmospheric  pressure  and  no  turbulence; 
j = 15  ma  = AU2 

where  typically  U - 10^  cm/sec. 

Then 

j = .015  amps  = A(10^)2 
A = 1.5(10”^)  (amp-sec2/cm2) 

and 

j = 1.5(10"10)U2  (31) 

Now,  equating  (30)  and  (31), 

nQ  = 90U2  (ions/cm3)  (32) 
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This  yields  a value  of  nQ  = 10-10  ions/cm  , which  is  a 
reasonable  value. 

Now,  plugging  this  value  of  nQ  back 

into  Eq.  (24), 

trec  ° 1/2(10‘6)(9/90U2)  and 

xrec  2 5(104)/U  (cm)  (33) 

4 

which  yields  for  U = 10  cm/sec, 

o c 

xrec  ^ cm 

It  can  be  seen  that  the  order  of  magnitude 
of  the  terns  involved  appears  to  approximately  agree  with 
the  observed  data.  This  tends  to  substantiate  the  basic 
assumptions  regarding  the  secondary  filaments  as  a recom- 
bination phenomenon.  Figure  (22)  shows  the  comparison  of 
Eq.  (33)  with  observed  data. 

Assuming  that  the  dependence  illustrated 
in  the  preceding  paragraphs  holds  for  the  "ordinary"  flow 
mode,  it  can  be  seen  that  this  idea  can  also  be  applied  to 
the  "virtual"  anode  concept  in  the  following  manner.  Since 
xrec  is  a measure  of  the  maximum  extent  of  the  space  charge, 
it  must  be  proportional  to  6 , the  anode  extension,  which 
is  also  a measure  of  space  charge  extent.  This  leads 
immediately  to  the  following  result: 
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*rec  * 5 * <1/U>C 

and 

vb  * Vbo(d- 6 ) = Vb0(d-C/U)  (34) 

so  as  U increases,  increases.  This  effect  has,  of 
course,  been  observed  as  noted  in  the  previous  section  at 
velocities  above  = 75m/sec.  As  a result  of  the  above,  it 
is  difficult  to  -determine  at  this  point  which  process, 
detachment  of  space  charge  from  the  anode  or  increasing 
recombination  rate,  or  a combination  of  the  two,  is 
responsible  for  the  upward  swing  of  the  V^/V^  vs  velocity 
curve  in  the  ordinary  flow  mode  (Fig.  23). 

(8)  Reverse  Flow  Convection  effects.  In  this 
orientation  (Fig.  9b),  with  no  flow,  the  static  anode  ex- 
tension exists  as  in  Figure  (17a) . As  velocity  is  increased, 
any  positive  ion  located  anywhere  off  the  pin  tip  stagnation 
streamline  will  experience  velocity  components  in  the  y 
direction  due  to  the  radial  component  of  the  anode  electric 
field  and  in  the  negative  and  positive  x directions  due  to 
convection  and  axial  electric  field  components.  The  total 
velocity  in  vector  form  will  be 
v = k+E+U 

where  the  positive  direction  is  from  anode  to  cathode. 


67 


- ..UBJJ 


Due  to  the  radially  directed  field  lines 
near  the  pin  tips,  the  resulting  space  charge  pattern  will 
tend  toward  the  pattern  of  flow  produced  with  a source 
flow  immersed  in  a uniform  flow  field.  The  result  is  that 
velocity  has  a direct  shortening  effect  on  the  anode  ex- 
tension; but  also  note  that  the  space  charge  tends  to  spread 
radially  which  reduces  the  charge  density  and  therefore  the 
field  strength.  The  reduction  of  6 and  spreading  tendency 
reduce  the  probability  of  streamer  formation  due  to  high 
charge  concentrations,  and  thus  reduce  the  current  flow. 
Reduction  in  current  means  less  charge,  and  so  the  glow 
collapses  (this  point  will  be  covered  in  more  detail  in  a 
later  Section) , which  quickly  reduces  the  space  charge 
distortion  to  the  point  that  the  breakdown  voltage  approaches 
the  upper  limit  as  <5  approaches  zero  (Fig.  24).  It  would 
be  logical  to  assume  that  most,  but  probably  not  all,  of 
the  space  charge  anode  extension  would  be  destroyed  or 
reduced  below  the  critical  level,  and  ^A^^max  = d/(d-  gQ) . 

It  can  be  seen  from  Figure  (14)  for  the 

reverse  flow  mode  that  this  was  the  case.  V./V,  ) for 

b bo  max 

a 3.9  cm  gap  is  predicted  by  the  present  model  as  1.16  and 
the  experimental  value  was  1.13  = v^^bo^max‘ 
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Werner  (Ref.  15)  reports  on  some  work 
with  a negative  point-plane  gap  with  velocity  from  the 
plane  (anode)  to  the  point  (cathode) . The  resulting  space 
charge  pattern  and  dynamic  processes  should  be  similar  to 
that  described  above,  only  with  the  negative  ions  playing 
the  part  of  the  positive  ions  in  the  former  case.  It 
might  be  expected  then  that  the  negative  point  would 
experience  trends  in  the  same  direction  as  a positive 
point  in  regard  to  breakdown  voltage  and  current  changes 
with  increasing  velocity.  It  turns  out  that  this  is  the 
case.  Breakdown  voltage  increased  with  increasing  velocity 
with  a maximum  V^/V^  = 1.08.  Calculating  V^/V^  )max  for 
this  case  from  Eq.  (20)  yields  ==  1.25,  so  the  voltage  ratio 
is  less  than  the  maximum.  Reference  (15)  also  reported  no 
distortion  in  the  glow  pattern  with  flow  while  being  ob- 
served under  a microscope  except  for  a reduction  of  size. 
This  would  be  a logical  result  of  the  process  described 
above  for  the  positive  point-to-plane  since  convection 
causes  glow  collapse,  and  an  equilibrium  glow  region  would 
be  established  which  was  much  smaller  than  the  no-flow 
glow. 

The  apparent  reduction  of  current  at  con- 


stant voltage  in  the  reverse  flow  mode  is  of  further 
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importance.  Reference  (16)  reported  no  appreciable  change 
in  breakdown  current  in  the  reverse  flow  mode  over  the 
ordinary  mode.  This  implies  a reduction  in  current  at  a 
given  voltage,  since  increased  appreciably;  while 
Reference  (13)  definitely  experienced  lower  currents  at  a 
given  voltage  with  flow  as  the  present  analysis  would 
expect.  This  fact  will  make  it  possible  to  determine  the 
correct  form  for  the  functional  dependence  of  the  dis- 
charge current  on  velocity  which  will  be  undertaken  in 
the  next  section. 

(9)  Analytic  Considerations  in  I = f (U) . While 

many  investigations  have  derived  expressions  for  current  vs. 

o 

velocity  of  the  form  I = U (Ref.  14) , it  can  be  shown  that 
a more  appropriate  expression  is  I * (exp)^. 

In  order  to  see  this,  it  is  necessary  to 
examine  some  facts  about  the  functional  observed  dependence. 
It  is  known  that 

I =>  f(U)  Condition  (1) 

and  within  practical  limits  of  velocity, 

1^0  Condition  (2) 

Also,  when  U = 0,  I = IQ  (the  no-flow  value).  Condition  (3) 
It  has  been  shown  also  that  when  U < 0 (reverse  flow)  the 
current  tends  to  decrease.  This  implies  that 
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dl/dU  ^ 0 at  U ■ 0 Condition  (4) 

and  from  References  (15,  16), 

dl/dU  - 0 at  U<0  Condition  (5) 

There  are  several  functions  which  could 
satisfy  conditions  (1  & 2),  but  two  of  the  simplest  are: 


(35) 


To  satisfy  condition  (3),  these  can  be  written: 

I-,  = AU2  + I 
1 o 

I2  = Io(exp)CU 


(36) 


Now, 


dIL/dU  = 2AU 


(a) 


CU  (37) 

dI2/dU  = IQC(exp)  (b) 

It  can  be  seen  that  Eq.  (37a)  does  not  satisfy  condition 
(5),  but  (37b)  does  provided  that  C is  a positive  number. 
Therefore,  the  equation  which  meets  all  the  appropriate 
conditions  can  be  written 

I = Io(exp)CU  (38) 

In  a later  section  it  will  be  shown  that, 
from  a physical  standpoint,  the  exponential  form  for 
velocity  dependence  evolves  quite  naturally,  although  in 
a slightly  different  form  than  Eq.  (38). 
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B.  GASDYNAMIC  EFFECTS  ON  DISCHARGE  CURRENT  (NON-UNIFORM 

FIELD) 

As  stated  in  the  Introduction,  a point- to-plane  dis- 
charge can  be  divided  into  distinct  regions  due  to  the 
divergence  of  the  applied  electric  field  distribution. 

These  are  (1)  an  intense  region  of  ionization  near  the 
point,  (2)  a space  charge  region  fed  from  the  ionization 
region  by  electrostatic  repulsion  of  positive  ions,  and 
(3)  the  quasi-neutral  region  in  which  the  applied  electric 
field  has  dropped  below  the  value  necessary  for  appreciable 
ionization  activity.  This  last  region  is  disturbed  more 
and  more,  however,  as  gap  breakdown  is  approached  by  the 
intrusion  of  pre-breakdown  streamers  whose  maximum  lengths 
are  proportional  to  applied  voltage  (Ref.  9).  The  actual 
length  of  any  one  streamer  is  a statistically  varying 
quantity  due  to  the  nature  of  the  dynamic  processes  in- 
volved. The  amount  of  current  associated  with  a specific 
streamer  is  proportional  to  its  length,  and  current  is 
also  proportional  to  streamer  frequency.  Since  both 
streamer  length  and  frequency  are  proportional  to  applied 
voltage,  the  current  must  be  nearly  proportional  to  the 
square  of  the  voltage. 

Several  empirical  expressions  have  been  used  to  des- 
cribe the  current-voltage  relationship  in  a point-to-plane 
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discharge  under  no-flow  conditions.  Reference  (7),  for 
example,  gives  the  following  empirical  formula  for  a posi- 
tive point -to -plane: 

I0-CV(V-M)  (39) 

where  M is  the  starting  potential  of  the  corona  current, 

and  C is  a function  of  the  geometry  and  gas  related 

o 

variables.  This  formula  gives  the  V dependence  previously 
mentioned,  and  has  been  found  to  satisfactorily  describe 
the  I-V  characteristic  under  no-flow  conditions. 

Combining  this  expression  with  the  current -velocity 
dependence  developed  in  the  preceding  section,  the 
equation  is  now 

I(U)  = I0(exp)CU 

where  IQ  is  given  by  Eq.  (39).  IQ  can  also  be  written  as 
IQ=CkV(V-Vc)  which  explicitly  reveals  the  density  dependence 
in  the  mobility  k=k0(  p/ p Q)  (Ref.  11) . This  expression 
is  also  too  limited,  however,  since  it  does  not  provide 
any  means  of  introducing  the  effects  of  turbulence  on  dis- 
charge current.  Since  turbulence  has  proven  to  be  a 
significant  parameter,  an  equation  to  describe  gasdynamic 
effects  on  a discharge  must  be  capable  of  reflecting 
turbulence  as  well  as  convection  and  density  change  effects. 
Such  an  expression  can  be  found  by  going  back  to  the  basic 
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physical  mechanism  of  current  production  in  a gaseous 
discharge,  the  electron  avalanche. 

In  the  special  case  of  a uniform  field  with  no  space 
charge,  the  current  (i)  produced  by  an  electron  avalanche 
can  be  expressed  as  a function  of  the  Townsend  first 
ionization  coefficient  a in  the  following  way  (Refs.  4,  7, 


8,  12): 


i = iQ(exp) 


This  equation  indicates  the  exponential  increase  in  current 
due  to  the  nature  of  the  electron  avalanche  process  and  is 
well  described  in  the  references  cited.  This  expression 
becomes  more  complicated  in  a non-uniform  field  where  a 
being  a function  of  the  electric  field,  changes  with  x. 

It  is  then  necessary  to  integrate  a over  the  distance  of 


avalanche  travel: 


i = iQ(exp) 


/. 

•'O 


U)dl 


For  purposes  of  the  present  discussion  an  effective 


can  be  defined  over  the  avalanche  travel  distance: 


aeff=(l/x)70  a(C)dt 


fo  a< 


where  x is  the  avalanche  travel  and  0=s  £ -x,  then 


i = iQ(exp) 


aeffx 
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The  subscript  notation  for  a will  be  dropped  for  conven- 
ience, but  it  must  be  remembered  that  in  a non-uniform 
field  with  space  charge,  a can  vary  in  the  space  of  a 
mean  free  path. 

A commonly  used  semi-empirical  expression  for  a is 
the  following  (Ref,  4) : 

a/p  = A(exp) ^ “BP/E)  (43) 

where  A and  B are  constants  of  the  specific  gas  under 
consideration.  For  air.  Reference  (4)  gives 
A = 15  ionizations/cm-torr 
B = 365  volts/cm-torr 
with  a reduced  field  range  of  validity  of 
E/P  = 100-800  volts/cm-torr. 

Of  course,  the  current  measured  in  a point -to -plane 
discharge  circuit  is  the  sum  of  many  small  avalanches  which 
form  primarily  near  the  point  in  the  high  f^eld  region; 
but  streamers,  which  propagate  far  into  the  low  field 
region  are  also  built  up  from  numerous  secondary  avalanches 
feeding  into  the  tip  of  an  advancing  streamer.  It  seems 
apparent  that  any  analysis  of  gasdynamic  effects  on  dis- 
charge current,  and  hence  power,  must  begin  at  this  level 
to  determine  how  the  gasdynamics  affect  the  electron 
avalanche. 
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. Pure  Convection  Effects  on  Current 

Consider  first  the  effects  of  pure  convection. 

The  model  used  in  discussing  breakdown  voltage  predicts 
that  extension  of  the  space  charge  region  will  initially 
occur  as  gas  velocity  is  increased.  This  anode  extension 
effectively  moves  the  boundary  of  the  high  field  region 
further  from  the  anode,  resulting  in  a longer  average 
electron  avalanche  travel  distance  (Fig.  25).  (Note  that 
in  discussing  breakdown  voltage,  the  apparent  anode  ex- 
tension was  thought  of  as  a reduction  of  effective  gap 
length,  and  in  discussing  current  effects,  it  is  an  in- 
crease in  avalanche  travel.)  Recalling  from  the  numerical 
example  computed  earlier  the  no -flow  anode  extension  60 
was  0.45  cm,  and  the  apparent  increase  above  this  due  to 
convection  was  0.15  cm.  An  avalanche  commencing  at  the 
beginning  of  the  critical  space  charge  region  with  con- 
vection would  only  travel  0.15  cm  further  than  in  the  no 
convection  case.  This  does  not  seem  like  much,  but  a 
simple  calculation  of  avalanche  current  from  Eq.  (40)  for 
both  cases  leads  to  surprising  results.  Using  an  average 
field  strength  of  Eavg=50,000  volts/cm,  P=760  torr: 

a = AP(exp)^  B^/^=44.4  ionizations /cm 
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and  since  (i/iQ)  - expa  X,  for  x = 0.45  cm  (no  flow); 

I/Vnf  " 4.77(109) 
and  for  x = 0.45  + 0.15  (with  flow): 

i/i0)F  = 3.7K1011) 

and  the  increase  of  current  for  this  single  avalanche  is 

ttH1-  4 80 

1/:lo'nf 

Although  the  numbers  used  above  are  only  order  of  magnitude 
correct,  the  obvious  conclusion  is  that  even  very  small 
distortions  of  E-field  structure  by  space  charge  convection 
can  cause  large  changes  in  current. 

It  is  worthy  of  note  at  this  point  that  the  expon- 
ential current-velocity  dependence  derived  in  the  previous 
section  from  purely  empirical  and  functional  considerations, 
turns  out  to  be  the  same  form  the  physical  model  would 
predict,  that  is  an  exponential  increase  due  to  increased 
avalanche  travel  in  the  high  electric  field  region. 

2.  Current -Vo ltage* Characteristics 

As  a further  test  of  the  space  charge  model  pro- 
posed, it  is  instructive  to  look  at  other  observed  I-V 
characteristics.  For  example,  the  standard,  no-flow  I-V 


curve  is  known  to  be  well  described  by  the  empirical 
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expression  of  the  form  of  Eq.  (39).  Now,  in  terms  of 
physical  processes,  consider  point-plane  electrodes  with 
increasing  voltage.  There  will  be  no  ionization  activity 
until  the  point  reaches  the  potential  required  for  ioniza- 
tion of  the  surrounding  gas  molecules,  which  gives  rise 
to  the  "inception"  voltage  (M  or  V£)  , at  which  the  corona 
current  starts  as  small  weak  pulses  caused  by  the  initial 
formation  of  the  electron  avalanches.  As  voltage  continues 
to  rise,  the  high  field  front  of  the  electric  field  moves 
out  into  the  gap,  causing  the  average  distance  at  which 
avalanches  can  form  to  move  with  it.  As  already  demon- 
strated, this  results  in  a rapid  rise  in  current,  which, 

from  a physical  standpoint,  is  better  described  by  an 

2 

exponential  function  rather  than  a U type  of  empirical 
fit. 

A simplified  representation  can  be  obtained  for 
the  total  current  by  observing  the  general  form  of  the 
result  of  a summation  of  several  individual  avalanche 
currents.  This  will  be  useful  in  later  developments  and 
for  a general  trend  analysis  to  further  test  the  validity 
of  the  proposed  model. 

Assume  that  the  total  number  of  avalanches  per 
unit  time  is  N.  Then  the  total  current  is 
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N 

I - Y i.  = i,  + i2  + ... 

J-l 

Since  most  activity  takes  place  in  the  same  small  region 
of  the  pin  tip,  let  the  following  average  quantities  be 
defined: 

N 

x = (1/N)  Y xi 

J-l 

N 

a = (1/N)  Y a i 

J-l  J 

N 

I - (1/N)  E iOJ 

J-l 

Summing  up  all  the  avalanche  contributions,  the 
total  current  is 

IQ-i1  + i2  + • • • = ° X + ^(exp)  “ X + . • . 

Io"ci(exp)  a X - C2 

where  C,  = i N and  C„  is  a function  of  the  electrode 
1 o 2 

geometry,  since  no  exponential  increase  in  current  can 
take  place  until  the  electric  field  strength  at  the  electrode 
surface  at  least  reaches  the  ionization  strength  of  the  gas 
surrounding  it.  Discounting  normal  saturation  currents 
due  to  free  electrons  which  are  always  present  in  the  air, 
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C2  can  be  evaluated  by  the  fact  that  I = 0 at  the  inception 


of  the  corona  current: 


!0  = 0 = C1(exp) 


C2  = C^(exp) 


at  corona  ignition 


but  at  the  point  when  ionizing  activities  first  commence, 


x = 0 , and 


C2  = C^(exp)  = = C 


and,  therefore, 


IQ  = C(expa  x -1) 


Equation  (44)  is  the  basic  physical  result  for  non- 
uniform  fields  under  no-flow  conditions  as  derived  from  the 
space  charge  model.  The  associated  constant  is  again  a 
function  of  gas  type,  density  and  electrode  geometry.  It 
is  interesting  to  note  that  this  is  nothing  more  than  the 
expression  for  current  between  parallel  plates  in  a uniform 
field,  written  to  accommodate  variable  geometry. 

Closer  inspection  of  the  constant  C = iQN  shows 
that  it  is  still  a function  of  streamer  length  L = f(V) 

O 

and  frequency  N = f = V,  so  can  be  rewritten  as  CV  and 
then  Eq.  (44)  can  be  written  as 


IQ  = C\T (exp  ax  -1) 
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or,  taking  one  V inside, 


a x 

I - CV (Vexp  - V) 

and  from  the  condition  that  1=0  and  x = 0 at  V = Vc , 
this  results  in 

a x 

Io  = CV(Vexp  -Vc) 

This  is  nearly  identical  with  the  previous  empirical  forms 
for  the  no-flow  condition,  except  that  now  the  exponential 
expression  introduced  permits  easy  interpretation  of  gas- 
dynamic  effects  on  the  discharge  current. 

The  x in  Eq . (44)  is  proportional  to  6 , so  it  will 
be  less  confusing  to  rewrite  the  current  equation  in  the 
following  form: 

IQ  = CV(VexpC  a 62vc)  (45) 

remembering  that  the  subscript  'o'  represents  no-flow 
conditions.  Equation  (44)  is  plotted  in  Figure  (26), 

It  is  now  easy  to  see  the  effects  of  increasing 
velocity  on  the  current,  since  in  the  ordinary  mode, 
increases  with  flow.  This  would  increase  the  exponent 
which  has  the  effect  of  shifting  the  I-V  characteristics 
up  and  to  the  left  as  in  Figures  (27  a&b) . These  curves 
are  typical  of  the  characteristic  shifts  which  have  been 
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observed  in  the  present  experimental  work.  The  general 
form  of  6 in  terms  of  the  other  variables  is 
<$  = <50  + UT  - f(u)  = c—  + UT  -f (u) , where  f(u)  is  the 

result  of  space  charge  "blowout"  and  recombination  effects. 

Looking  again  at  the  reverse  flow  case  discussed 
earlier  with  the  aid  of  Eq.  (45),  it  is  now  easy  to  see 
how  a low  current  equilibrium  situation  is  established. 

As  soon  as  there  is  even  a slight  distortion  (shortening) 
of  the  space  charge  region,  the  current  decreases  along  an 
exponential  curve,  which  results  in  less  charge  production, 
further  decrease  in  current,  etc.,  until  equilibrium  is 
again  established,  when  the  space  charge  region  shrinks 
far  enough  into  the  high  field  region  of  the  anode  that 
it  is  no  longer  affected  by  convection.  Figure  (28) 
summarizes  the  effects  of  velocity  on  6 and  the  breakdown 
voltage  characteristic  curves. 

3.  Diffuser/Density  Effects  on  Current 

The  effects  of  reduced  density  on  the  current  are 
also  easily  deduced  from  Eq.  (45),  since  as  shown  earlier, 

6 increases  with  a density  reduction,  but  also  a ( p) 
increases  with  decreasing  density  over  the  ranges  of 
reduced  pressure  discussed  so  far.  (There  is  an  optimum 
1 n Lzation  efficiency,  n a/E,  which  occurs  for  air  at  an 


82 


E/P  * 365  volts/cm-torr  = B,  which  will  be  discussed  below 
(Ref.  12)).  The  result  of  this  dual  exponential  effect  is 
a rapid  rise  in  current  with  decreasing  density,  until  the 

f 

reduced  pressure  in  the  ionization  region  drops  below  the 
optimum.  Beyond  that  point,  the  effects  of  density  reduc- 
tion on  the  current  could  actually  reverse  if  the  ionization 
efficiency  drops  too  far  (Fig.  29). 

To  arrive  at  an  estimate  for  optimum  pressure  for 
best  ionization  efficiency,  it  is  possible  to  write 


E/P)  t=  365  v/cm-torr  = B 

§ 

and  using  E = 200,000  volts/cm  (value  at  pin  for  a 3 cm 


gap,  hyperboloidal  point-to-plane) , 

■ 

popt=  E /B  = 200,000/365  - 547  torr 


The  value  computed  from  the  data  of  Reference  (18) 
in  Appendix  (b)  for  a 3 cm  gap  is  = 480  torr,  which  is  not 
too  far  off,  considering  the  approximation  inherent  in  the 
point  model  used.  (For  example,  if  the  point  radius  were 
ar jumed  to  equal  the  wire  radius  of  = 0.08  cm,  then 
P0pt=375  torr) . Since  the  discharge  in  the  present 
investigation  was  operated  at  a pressure  well  above  the 
optimum  value  for  most  efficient  ionization,  it  would  be 
expected  that  any  density  reduction  produced  in  the  test 
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section  by  the  addition  of  a diffuser  would  enhance  ioniza- 
tion significantly,  which  it  does  (Fig.  30). 

4.  Turbulence  Effects 

While  the  study  of  turbulence  is,  in  itself,  a 
complicated  field,  the  basic  interactive  effects  of  turbu- 
lence on  an  electric  discharge  can  be  qualitatively  des- 
cribed with  only  the  most  basic  characteristics  of  a 
turbulent  flow  field;  those  of  fluctuating,  time-dependent 
variations  in  gas  density  and  velocity. 

The  turbulence  generated  in  the  discharge  section 

o 

in  the  present  work  ranged  in  frequency  from  = 100-20,000  Hz 
(Ref.  3).  This  frequency  range  yields  characteristic  times 
t = (1/f)  of  the  order  of  2(10"^)-5(10  ~>)  sec.  These 
times  are  long  compared  with  the  characteristic  times 
associated  with  charged  particle  kinetic  processes  (Table  I) . 
The  result  is  that  the  charged  particle  properties  adjust 
in  a quasi-steady  fashion  to  fluctuations  in  gas  properties. 
Thus,  charged  particle  kinetics  and  energy  transfer  pro- 
cesses couple  to  fluctuations  in  the  neutral  gas"  primarily 
through  changes  in  gas  temperature  and  density  (Ref.  19). 

For  example,  consider  a point  in  a turbulent  dis- 
charge region;  gas  density  at  this  location  will  vary 
periodically  around  a mean  value.  Although  the  time 
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average  of  the  fluctuations  at  the  point  is  zero,  and  the 
point  experiences  as  much  time  at  elevated  pressures  as 
it  does  at  reduced  pressures,  the  net  effect  of  these 
fluctuations  on  discharge  current  are  very  significant. 

This  is  due  to  the  exponential  dependence  of  the  local 
value  of  the  ionization  coefficient  on  the  local  reduced 
pressure  (E/P).  As  long  as  the  discharge  (E/P)>  B, 
pressure  fluctuations  will  significantly  increase  the  aver- 
age effective  ionization  coefficient. 

EXAMPLE:  P = 760  torr  = mean  pressure 

p = 50  torr  = fluctuating  component 
a = AP(exp) (-BP/E) 

and  using  a simple  expression  for  current  for  illustration: 
I = I0(exp)aX 

Computing  ionization  coefficients  for  E = 50,000  v/cm: 


“760 

= 44.4 

Mean 

“810 

= 32.9 

High  P fluctuation 

“710 

= 59.8 

Low  P fluctuation 

Current  computations  for  I = 1,  x = 0.5  cm: 

I = (exp)^*^(*5)=  4.38(10^)  Mean 

I = (exp)32'9(,5)=  1.39(107)  High  P 

I = (exp)59,8(*5)=  9.67(1012)  Low  P 

Iavg=  tthi  + Ilo>/2  " 4.83(1012) 
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The  amplification  of  the  current  over  the  mean  due  to  a 
fluctuation  in  pressure  is  then: 

WImean-4-83<1°12)/4-38<1°9)  * 1103 

Obviously,  the  effects  of  turbulent  fluctuations  in  pressure 
or  density  can  cause  substantial  increases  in  current  due 
to  the  quasi-steady  coupling  of  electron  kinetics  to  temporal 
variations  of  the  neutral  gas  properties. 

Local  increases  in  electron/current  density,  however, 

. 

can  have  a great  destabilizing  effect  on  the  discharge, 
since  increased  local  current  density  causes  increased 
joule  heating  ( j * E)  of  the  gas.  This  heating  results  in 
local  density  reduction,  increased  conductivity  and  further 
increases  in  current  density  which  may  result  in  a runaway 
condition  leading  to  complete  breakdown  of  the  gap  (Ref.  20) . 

Fortunately,  it  is  the  basic  nature  of  turbulence 
that  the  strongest  density  fluctuations  are  accompanied  by 
the  highest  fluctuating  velocity  components  (Ref.  19) . 

These  high  velocity  components  tend  to  disperse  concentra- 
tions of  space  charge,  dissipate  local  hot  spots  and  pro- 
mote stability.  The  role  of  turbulence,  then,  is  two-fold: 
on  the  one  hand,  it  enhances  ionization,  while  on  the  other 
hand,  it  jointly  acts  to  prevent  an  unstable,  runaway 
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condition  from  developing  due  to  formation  of  local  hot 
spots  in  the  gas. 

Before  leaving  the  subject  of  turbulence,  it  is 
worthwhile  to  do  a general  frequency  analysis  to  determine 
how  a given  turbulent  frequency  spectrum  will  affect  a 
discharge.  It  was  pointed  out  in  Reference  (16)  that 
turbulence  spectra  which  were  heavily  weighted  in  the  low 
frequency  regions  gave  the  best  improvement  in  power  enhance- 
ment to  the  discharge. 

That  this  should  be  so  can  be  seen  by  analyzing  the 
effects  of  turbulence  more  explicitly  on  the  ionization 
coefficient  a . In  a turbulent  field,  a = cx(t),  and  if 
f(co  t)  is  a periodic  fluctuation  of  gas  density,  then  the 
time  average  of  a(t)  over  one  cycle  will  be: 

a(t) dt 

but  since  co  is  an  independent  parameter  of  the  problem, 
it  may  be  included  as  follows: 

(nit)d((ot)  = — (46) 

CO 

This  shows  that  a is  inversely  proportional  to  the  turbu- 
lent frequency  components.  It  should  be  possible  to 
develop  a weighting  factor  for  a given  turbulence  spectrum 


a = 1/co  T J 
Jo 
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by  svimming  the  contributions  of  the  various  frequency 
components  to  the  ionization  coefficient.  It  must  be 
realized,  however,  if  an  attempt  is  made  to  quantify  the 
effects  of  a given  turbulence  spectrum,  that  inhomogenieties 
in  the  turbulence  pattern  can  lead  to  premature  breakdown 
through  a relatively  Inactive  "dead"  region  of  the  gap 
where  instabilities  may  develop  uninhibited. 

The  following  is  a summary  of  the  gasdynamic  effects 
on  discharge  current  observed  in  the  present  experimental 
work  and  interpreted  in  terms  of  the  present  model: 

a.  No -Flow 

I = CV  (VexpCia<S°-V-) 
o \ c 

Current  commences  at  an  inception  voltage  Vc 
and  rises  on  a modified  exponential  curve  until  breakdown 
occurs  (Fig.  26) . 

b.  Flow  (ordinary  mode) 

I(U)  = CV(VexpCa5  (U)-Vc) 

The  anode  extension  6 increases  due  to  con- 
vection of  space  charge,  shifting  the  characteristic  curve 
up  and  left,  while  a also  increases  due  to  high  space 
charge  in  the  gap  (Fig.  27b). 
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c.  Reduced  Density  (diffuser  on) 

6 and  a increase  due  to  reduced  density,  and 
6 also  increases  due  to  the  convection  effect.  The  result 
is  a large  shift  of  characteristics  up  and  left.  The  re- 
duced density  also  reduces  the  breakdown  voltage  signifi- 
cantly, and  may  result  in  a reduction  of  maximum  power 
(IbVb)  (Fig.  30). 

d.  Turbulence  (diffuser  off) 

Substantial  improvements  in  a can  be  realized 
with  the  proper  type  of  turbulence,  and  stability  is  con- 
currently enhanced.  Mean  velocity  also  causes  6 exten- 
sion due  to  convection,  and  a large  increase  in  a 6 
results,  with  the  accompanying  shift  up  and  left  of  the 
I-V  curve  (Fig.  27a). 

e.  Turbulence  and  Density  Reduction  (diffuser  on) 

The  combined  effects  of  (c)  and  (d)  above 

result  in  maximum  I-V  shift  and  high  power  capacity  due  to 
stability  enhancement  by  turbulence. 

Figures  (19,  20)  are  a sequence  of  pictures  of 
the  discharge  section  and  the  discharge  operating  under 
some  of  the  conditions  described  above,  and  Figure  (18)  is 
a schematic  of  the  pin  tip  regions  as  they  appeared  during 
visual  observation.  One  additional  feature  of  turbulence 
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which  is  evident  from  these  figures  is  that  in  addition  to 
homogenizing  the  discharge  region  the  turbulence  creates  a 
much  larger  ionization  region  at  the  anode,  which  is  the 
major  source  of  the  additional  current. 

5 . Maximum  Current/Discharge  Stability 

Figure  (31)  shows  how  the  breakdown  current  varied 
in  the  present  investigations  as  a function  of  mass  flow 
rate  through  the  discharge  region  for  the  various  con- 
figurations already  mentioned. 

The  curve  for  pure  convection  (no  diffuser } no 
turbulence  grid)  shows  increasing  stability  with  increasing 
mass  flow  rate.  This  stability  increase  is  probably  due, 
at  least  in  part,  to  convection  of  heat  out  of  the  high 
current  density  regions  near  the  pin  tips  where  growing 
thermal  instabilities  may  develop  under  no-flow  conditions. 
The  data  indicate  that  a limit  is  apparently  being  approached 
in  this  mode.  This  would  be  expected  in  a relatively  smooth 
laminar  type  of  flow,  since  the  joule  heat  developed  at  a 
pin  tip  would  flow  toward  the  cathode  with  relatively  slow 
radial  diffusion,  resulting  in  a column  of  heated  gas  be- 
tween anode  and  cathode,  increasing  the  conductivity  and 
breakdown  probability  of  the  gap. 
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The  addition  of  the  diffuser  raises  the  curve, 
thus  indicating  a stabilizing  effect  associated  with  re- 
duced density.  This  stabilization  can  most  probably  be 
traced  to  the  increased  diffusion  rates  of  charged  parti- 
cles, reducing  high  density  charge  build-up,  and  increased 
thermal  conductivity,  which  also  tends  to  reduce  local 
heating.  The  problem  of  laminar  flow  is  still  present, 
however,  and  it  can  be  seen  that,  once  again,  a current 
or  power  limit  is  being  approached  as  mass  flow  rate 
increases . 

The  addition  of  the  turbulence  generating  screen 
in  either  of  the  above  configurations  demonstrates  a very 
marked  increase  in  stability,  with  no  apparent  saturation 
limit  being  approached  within  the  limits  of  the  power 
supply  capabilities  used  in  the  present  investigations. 

It  thus  appears  that  the  fluctuating,  large  scale,  high 
intensity  velocity  components  introduced  are  very  effective 
in  enhancing  discharge  stability.  This  enhancement  is 
presumed  to  act  through  intense  mixing  of  ionic  space 
charges  with  surrounding  areas,  and  rapid  dissipation  of 
heat  from  local  regions  of  the  discharge  gap.  The  result 
is  a smoothing  out  of  the  electric  field  distortion  caused 
by  space  charge,  and  a rapid  redistribution  of  heat 
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throughout  the  entire  discharge  volume,  whereas  in  the 
laminar  flow  cases,  the  heat  was  confined  to  relatively 
small  columnar  volumes  between  anode  and  cathode. 

The  most  effective  screen  design  turned  out  to  be 
a redesigned  "number  9"  screen  of  Reference  (6) . Reference 
(3)  gives  detailed  engineering  data  on  all  equipment  used 
in  the  present  work.  One  of  the  significant  features  of 
the  redesign  was  the  incorporation  of  large  boundary  layer 
bleed  areas  around  the  discharge  volume  (Fig.  2).  Flow 
velocities  measured  through  this  boundary  layer  region 
were  2-3  times  the  average  discharge  velocity  with  the 
turbulent  screen  in  place.  It  is  felt  that  this  high 
velocity  region  serves  as  a heat  sink  for  the  discharge 
and  greatly  enhances  the  thermal  stability  by  very  rapid 
removal  of  waste  heat  generated  in  the  discharge  region. 

Reference  (3)  should  also  be  consulted  for  an 
interesting  approach  to  power  input  optimization  for  the 
type  of  flow  device  used  in  this  work. 
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V.  CONCLUSIONS 


It  appears  that  the  space  charge  convection  model 
presented  here  is  a valid  approach  to  electric  discharge 
and  gasdynamic  interaction  analysis  in  highly  non-uniform 
fields.  The  model  has  also  been  shown  to  hold  as  an 
approximation  to  discharge  breakdown  properties  in  a 
uniform  field,  and  so  may  not  be  limited  to  non-uniform 
geometries  where  it  appears  to  describe  many  observable 
effects  concerning  gasdynamic  interactions. 

Gasdynamic  interactions  promote  stability  of  a dis- 
charge in  several  ways,  through  pure  convection,  density 
reduction,  and  turbulence.  It  appears  that  the  basic 
mechanism  of  all  three  of  these  is  through  thermal  stabi- 
lization by  heat  removal  from  potentially  developing  "hot 
spots"  in  the  gas. 

High  intensity,  low  frequency  turbulence  is  best  for 
stabilization.  This  fact  may  permit  quantitative  analysis 
of  individual  turbulence  spectra  for  use  in  optimization 
studies. 

The  addition  of  a diffuser  and  turbulence  generating 
screens  results  in  much  higher  power  input  at  lower  mass 
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flow  rates,  and  therefore  is  an  important  economic 
consideration  in  the  construction  of  larger  devices  of 
this  type. 
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VI.  RECOMMENDATIONS 


r 

I 

I 

( 

A.  The  major  recommendation  resulting  from  this  study 
is  that,  initially,  more  comprehensive  testing  should  be 
done  on  all  aspects  of  the  gasdynamic/f low  interactions 
reported  in  order  to  establish  accurate  empirical  inputs 
required  for  optimization  studies.  Concurrently,  attempts 
should  be  made  to  derive  accurate  relations  among  the 
parameters  discussed  using  the  Poisson  Equation  and 
appropriate  flux  equations. 

B.  The  stability  problem  is  probably  the  most  critical 
and  complex,  but  since  it  is  the  key  limiting  performance 
parameter  for  high  pressure  discharges,  an  accurate  means 
of  quantification  is  necessary  which  properly  reflects  the 
influence  of  gas  dynamics.  It  is  hoped  that  the  relatively 
simple  model  presented  here  will  aid  in  visualization  of 
the  stability  problem  and  the  gasdynamic  factors  which 
affect  it.  In  regard  to  discharge  stability,  it  is 
specifically  recommended  that  the  following  tests  be 
performed: 

1.  Current  distribution  studies  should  be  performed 
for  correlation  with  temperature  profiles  and  breakdown 
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location  statistics,  and  turbulence  and  velocity  profile 
measurements . 

2.  Oscillographic  studies  of  current  and  voltage  wave- 
form characteristics  should  be  attempted  in  the  various 
configurations  tested,  and  at  various  gap  lengths  to  deter- 
mine the  exact  type  of  breakdown  which  occurs . 

3.  Small  diameter  wire  should  be  tried  as  the  anode, 
with  the  wires  mounted  parallel  to  the  cathode.  Wire 
electrodes  resulted  in  a much  more  uniform  discharge  in 
recent  MHD  laser  cavity  experiments  (Ref.  22). 

C.  Since  the  present  work  was  limited  by  power  supply 
capabilities,  a larger  power  supply  with  much  greater 
power  output  is  necessary  to  extend  the  investigation  of 
optimum  gap  length,  and  maximum  power  input  limits. 

D.  An  attempt  should  be  made  to  determine  more  specif- 
ically the  mechanism  responsible  for  the  increase  in  at 
high  velocities  in  the  ordinary  flow  mode,  since  small 

increases  in  V,  at  that  point  in  the  curve  result  in  a 
b 

large  increase  in  1^  due  to  the  slope  of  the  I-V 
characteristic . 

E.  Extension  of  this  work  to  a lasing  mixture  should 
be  attempted  when  a better  power  supply  is  available.  It 
is  believed  that  a lasing  mixture  may  prove  more  thermally 
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stable,  since  most  of  the  input  energy  goes  into  the 
vibrational  molecular  states  instead  of  translational 
gas  heating,  and  therefore  even  more  power  can  be  stably 
coupled  into  the  gas. 
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Figure  6.  Electric  Field  and  Voltage  Distribution  between 
Parallel  Plates.  Subscripts:  (1)  No  Space  Charge 

(2)  Moderate  Space  Charge  (3)  High  Space  Charge 
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★ Bandel  4 cm  Gap 

O Bandel  8 cm  Gap 


Figure  8.  Predicted  Model  Breakdown  Voltage  Ratio 
Density  Dependence  vs.  Actual  Data. 
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Figure  11.  Voltage  Ratio  vs.  Velocity  in  Referse 
and  Cross -Flow  Modes. 


Figure  12.  Electric  Field  Distribution  in  a Point-Plane 
Gap.  Curves:  (1)  Point  Field  Only  (2)  Space  Charge 

Only  (3)  Total  Field 
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CORRECTED  DATA 


GAP  LENGTH  (cm) 

Figure  13.  Breakdown  and  Corona  Initiation  Voltage  vs.  Gap  Length 
for  No-Flow  at  Atmospheric  Pressure. 
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Figure  15.  Calculated  Values  of  the  Slip  Parameter  at 
the  Field  Point  E = 20,000  volts/cm. 
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Figure  16.  Calculated  Values  of  Slip  Parameter  at  Anode 
Surface.  Note  that  Sp  = f(d) . 
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( a)  NO  FLOW 


( b)  INTERMEDIATE  FLOW 


60  to  90  m/sec 


(C  ) HIGH  FLOW 


over  75  m/sec 
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Figure  17.  Various  Stages  of  Anode  Extension  and  "Blowout1 


( c)  FLOW,  WITH  TURBULENCE 


Figure  18.  Visual  Character  of  Various  Glow  Modes 
without  Diffuser 


Figure  19.  Photograph  of  Discharge.  No  Flow 

(10  sec  exposure).  Pins  on  Left.  Bright  region 
on  bottom  is  an  arc  which  took  place  during 
timed  exposure.  Note  small  spherical  tip  glow 
region. 
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Flow.  No  Turbulence 


Flow.  With  Turbulence 


Figure  20.  Pins  on  Left.  Ordinary  flow  from  left  to 
right  (1  sec  exposures).  Note  elongation  of  glow 
region  in  (a)  due  to  convection  (u  - 150  m/sec) , and 
very  homogeneous  spreading  of  glow  in  (b)  due  to  high 
turbulence. 
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Schematic  Representation  of  Cross-Flow  Current  Distribution  (Observations). 


Ordinary  Flow  Mode. 
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1 0 ) Electric  Field 

(1)  No  Space  Charge 

(2)  With  Space  Charge,  No  Flow 

(3)  With  Space  Charge,  and  Flow 


(b)  Schematic  of  Flow  Effect  on  Space  Charge 


Figure  25.  Convection  Effects  on  Average  Electron 
Avalanche  Travel  Distance 
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Figure  26.  I-V  Characteristic  Velocity  Dependence  from 


Figure  27.  Convection  and  Velocity  Effects  on  I-V 
Characteristics.  (a)  3 cm  gap,  with  Turbulence 
(b)  2.9  cm  gap,  no  Turbulence. 
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APPENDIX  A 


VARIABLE  VOLTAGE  RATE  TESTS 

As  stated  in  the  Introduction,  this  series  of  tests  was 
undertaken  in  an  attempt  to  determine  whether  any  type  of 
resonance  condition  might  exist  between  the  rates  of  voltage 
increase  used  and  any  of  the  dynamics  of  the  discharge  or 
gas  parameters.  The  fact  that  the  characteristic  times  for 
current  pulsing,  rates  of  voltage  increase  and  turbulence 
all  overlap  in  the  ranges  tested  was  the  primary  motivating 
factor. 

The  breakdown  voltage  was  found  to  decrease  in  a linear 
fashion  with  increasing  voltage  rate  ( Y ) in  the  ordinary 
and  cross-flow  modes  (Fig.  A-l) . At  higher  velocities  in 
the  cross-flow  mode,  the  breakdown  voltage  increased  at  the 
very  low  values  of  Y as  would  be  expected  from  the  space- 
charge  -convection  model.  The  slope  of  the  curves  remains 
basically  unchanged,  however. 

In  the  ordinary  flow  mode,  the  effect  of  the  flow 

definitely  increased  the  slope  of  the  V curve,  but  V still 

b b 

decreased  linearly  with  increasing  voltage  rate. 

Under  no-flow  conditions,  using  the  Townsend  type  of 
formula  for  the  relation  between  current  and  voltage,  it 
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can  be  shown  that  the  slope  of  these  vs.  y character- 
istics should  be  negative.  Assume  is  a linear  function 
of  y , and 

V»>  -vbo  + AY 

where  the  'o'  indicates  the  value  at  f = 0.  The  Townsend 
current  equation  for  breakdown  is 

!b  - cVVb  ' V 

and  breakdown  power  is  Pb  = I^V^. 

pb  - CTb<W  * C<Vbo+A  )3-C<Vb0+A  >2yc 

and 

dPb/d  = 3C(Vbo+A  )2A-2C(VbQ+A  )kVQ 

Now,  assuming  P^  has  a maximum  or  minimum,  letting  the 
above  = 0: 

/ • = (2V_/3  - V,  )/A 

max/mm  v c'  bo 

Since  is  always  positive  and  V 2W^/3,  then  A,  the 

slope  of  the  Vb  vs.  curve,  must  be  negative  as  observed. 

The  reason  for  this  decrease  in  V,  and  I,  is  not  known 

b b 

and  was  not  pursued  due  to  the  uninteresting  trends  observed. 
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Figure  A-l.  Variable  Rate  of  Voltage  Increase  Effect  on  Breakdown  Voltage  for 
Cross  and  Ordinary  Flow  Modes. 


APPENDIX  B 


NO -FLOW  EMPIRICAL  GAP  LENGTH/ DENSITY 
OPTIMIZATION  AND  COMPARISON  WITH  OBSERVED  DATA 

Prior  to  commencement  of  the  testing  phase  of  this  work, 
data  was  gathered  from  a number  of  sources  on  the  break- 
down properties  of  a point-to-plane  gap.  Empirical  equa- 
tions for  V^(d,  ) and  I^d,  ) were  developed  and  used  in 

a computer  optimization  study  for  this  electrode  configuration. 

The  results  of  this  study  indicated  that,  under  no-flow 
conditions,  maximum  power/cm  input  could  be  obtained  at 
d = 4 cm  and  a density  ratio  of  0.65  = ( / ) . This 

corresponds  to  a pressure  of  = 420  torr  (Fig.  B-l) . Figures 
B-2  & B-3  show  the  computer  generated  and  1^  character- 
istics plotted  against  density  ratio.  These  figures  indi- 
cate that  the  4 cm  gap  was  best  at  atmospheric  pressure 
too,  while  data  from  the  present  experiments  yielded  an 
optimum  gap  of  = 1.6  cm  under  no-flow  conditions  at  atmos- 
pheric pressure  (Fig.  B-4)  for  a multiple  pin  electrode. 

The  equations  used  in  the  optimization  study  are  listed 
below  for  reference: 

Breakdown  Voltage: 

V,  ( ,d)  = 51(  / )1,16(1.8-  / )(l-expT)  + B 

D O O 
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where 


T = d1,5/3 

B = 20(d/do-l)(  / )'65(1-  / ) 

d = 4 cm 
o 

Breakdown  Current: 

V ’d’V  * '09c((  1 )Eyb  ■ A) 

A = 590(d/do-l)(  / -.04) * 62(1- 2-  / ) 

E = 1.05  + ( / - .05) * 33 

C = 1 for  d 4 cm 

C = .lld^*^  for  d 4 cm 

A series  of  tests  in  the  ordinary  flow  mode  was  con- 
ducted using  the  turbulence  screen  number  "9"  (Fig.  2)  of 
Reference  (6)  at  various  velocities  and  gap  lengths  in  an 
effort  to  determine  the  optimum  gap  length  at  atmospheric 
pressure  in  the  present  apparatus.  Figures  B-5  through  B-8 
show  the  results  of  these  tests. 

In  Figure  B-5  the  typical  V vs.  velocity  trend  is  evi- 
dent, showing  an  initial  decrease  of  and  a later  increase, 
except  for  the  short  gap-high  current  case  of  the  1.6  cm 
gap  which  was  explained  earlier. 

Figure  B-6  is  interesting  in  that  it  shows  1^  reaching 
a maximum  value  of  = 17  ma  for  the  three  best  gaps.  This 


133 


k 


may  be  due  to  the  fact  that,  with  no  boundary  slots  as  in 
later  tests,  the  blockage  of  this  screen  is  = 607o,  and  maxi- 
mum velocity  is  limited  to  = 120  m/sec  through  the  discharge 
section.  This  may  result  in  the  section  becoming  thermally 
unstable  due  to  inefficient  heat  removal. 

The  maximum  specific  power  is  shown  in  Figure  B-7.  Here 
it  can  be  seen  that  the  short  gap  permits  highest  specific 
power.  This  would  be  expected,  since  the  shorter  the  gap, 
the  lower  the  maximum  temperature  rise  in  the  gas  at  a 
given  velocity  (Ref.  5) . 

Figure  B-8  shows  the  breakdown  current  variation  with 
gap  length  for  various  constant  velocity  profiles  from  the 
data  of  the  previous  figures.  It  is  interesting  that  at 
high  velocity  and  high  turbulence,  the  4 cm  gap  current 
increases  faster  than  at  the  other  gap  lengths. 

It  must  also  be  noted  that  it  was  later  discovered  that 
the  power  supply  used  in  these  tests  was  limited  to  18-19  ma 
by  internal  problems,  so  part  of  the  apparent  17  ma  limit 
in  the  discharge  illustrated  in  Figure  B-6  may  be  attributable 
to  a power  supply  limitation  rather  than  a discharge  stability 
limit. 
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Gao  Length 


Figure  B-2.  Breakdown  Voltage  vs.  Density  Ratio  for  Single  Point-Plane 


Bandel's  Data 


Breakdown  Current  vs.  Density  Ratio  for  Single  Point-Plane 


Breakdown  Voltage  vs.  Velocity  for  Multi-Pin  Configuration  with 
Turbulence  Screen  #9.  (Various  Gaps) 


Breakdown  Current  vs.  Velocity  for  Multi-Pin  Configuration 
Turbulence  Screen  #9 


Figure  B-7.  Breakdown  Power/cm  vs.  Velocity  for  Multi-Pin  Configuration 
with  Turbulence  Screen  #9 
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